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fluctuations  in  the  horizontal  component  of  the  earth's 
magnetic  field  were  measured  at  a  land  site  (La  Mesa)  and 
on  the  floor  of  Monterey  Bay  at  a  depth  of  70  meters.  The 
measurements  span  a  19-hour  period  (August  17/18,  1982). 
Assuming  a  circular  polarization  basis  set,  analysis  was 
conducted  to  determine  the  coherence  between  the  two  sites 
in  the  frequency  range  .05-10  Hz.  Preliminary  analysis 
indicates  a  high  degree  of  coherence  from  .05 -.4  Hz.  Also, 
a  study  of  the  polarization  of  the  micropulsation  field  at 
each  site  was  conducted  in  terms  of  Stokes  parameters. 
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I.  INTRODUCTION 


This  thesis  research  is  part  of  an  ongoing  effort  being 
conducted  by  the  Naval  Postgraduate  School  to  obtain  a  long¬ 
term  data  base  of  ELF  electromagnetic  noise  on  the  ocean 
floor.  The  result  of  this  project  should  be  a  better  under¬ 
standing  and  interpretation  of  that  noise  and  its  sources. 

Magnetic  noise  is  important  to  the  geophysicist  and  to 
the  Navy.  Magnetic  techniques  are  applied  in  exploration 
for  minerals,  the  prediction  of  communications  propagation 
and  to  the  modeling  of  the  sphere  on  which  we  live.  Geo¬ 
magnetic  studies  are  of  interest  to  the  Navy  for  their 
application  to  submarine  detection,  submarine  communications 
and  to  mine  warfare. 

The  objective  of  this  thesis  within  the  broader  context 
of  the  overall  project  was  to  design  and  validate  a  computer 
code  that  would  facilitate  data  manipulation  and  analysis. 

To  that  end,  the  present  computer  code  calculates  for  two 
sites  (one  on  the  ocean  floor,  the  other  on  land),  the 
individual  site's  power  spectral  densities  and  Stokes 
parameters  0,1,  2  and  3  for  two  of  the  individual  site's 
orthogonal  components  in  the  horizontal  plane.  The  code 
also  computes  the  coherence  of  these  two  orthogonal  compo¬ 
nents  and  their  phase  difference.  Lastly,  the  program 
computes  the  right  and  left  circular  polarization  power 
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spectral  densities  of  the  sites  and  the  coherence  between 
the  seafloor  and  land-based  sites. 

As  discussed,  the  present  code  applies  only  to  two 
sites  and  two  orthogonal  geomagnetic  components  at  each 
site.  However,  extension  to  several  sites  and/or  components 
should  be  an  easy  undertaking  due  to  the  style  of  the  present 
program. 
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II.  BACKGROUND 


A.  THE  MAIN  GEOMAGNETIC  FIELD 

Historical  observation  and  modern  measurement  demonstrate 
that  the  principal  source  of  the  geomagnetic  field  (approxi¬ 
mately  50,000  nanoteslas)  is  beneath  the  surface  of  the 
earth.  An  easy-to-picture  causal  interpretation  of  this 
phenomenon  is  painted  by  imagining  the  earth  with  a  dipole 
magnet  near  its  center  as  pictured  in  Figure  2.1.  The 
magnet  axis  is  tilted  at  an  angle  of  11.5  degrees  from  the 
earth’s  rotation  axis.  This  model  provides  a  prediction 
accuracy  to  a  few  percent  of  the  measured  geomagnetic  field. 
Figures  accurate  to  less  than  one  percent  of  the  earth's 
measured  field  can  be  obtained  using  more  sophisticated 
models  [Ref.  1]. 

The  generally  accepted  physical  process  for  the  genera¬ 
tion  of  the  geomagnetic  field  is  that  of  the  self-sustaining 
dynamo  [Ref.  2],  Figure  2.2  presents  a  simple  picture  of 
such  a  dynamo.  The  physical  requirements  of  this  theory 
are  mostly  met  in  the  interior  of  the  earth.  For  example, 
the  earth's  center  is  a  large  conducting  system.  Addition¬ 
ally,  we  can  assume  that  the  earth's  molten  core  is  of  a 
differential  rotational  nature  with  turbulent  upwelling. 

Energy  sources  for  such  a  dynamo  are  most  likely  avail¬ 
able  in  the  earth's  core  due  to  a  number  of  processes.  This 
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Figure  2.1 
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Eccentric  Dipole  Model  of  Geomagnetic  Field 
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Figure  2.2  Simple  Disk  Dynamo 

all  boils  down  to  mean  that  the  motion  of  a  conductor,  such 
as  the  earth’s  molten  core,  in  a  magnetic  field  produces  a 
current  which  in  turn  produces  a  magnetic  field  in  support 
of  the  original  magnetic  field  [Ref.  3]. 

A  simplified  picture  of  the  dynamo  theory  originates 
by  assuming  that  in  the  distant  past  a  magnetic  "seed’'  field 
was  captured  by  the  earth’s  conductive  core.  The  differen¬ 
tial  rotation  of  this  core  ”winds  up”  the  field  into  an 
intense  azimuthal  field,  such  as  would  happen  to  a  twisted 
rubber  band.  Turbulent  upwelling  of  the  core  tends  to 
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carry  this  azimuthal  field  outward,  converting  it  to  a 
toroidal  helical  field  within  an  associated  azimuthal 
circulating  current.  This  current  in  turn  generates  the 
earth's  dipolar  field. 

An  additional  advantage  of  the  self-sustaining  dynamo 
theory  lies  in  its  ability  to  explain  the  well -documented 
reversals  of  the  geomagnetic  field.  It  has  been  postulated 
that  two  or  more  coupled  "dynamos"  as  shown  in  Figure  2.3 
can  cause  reversals  of  the  main  geomagnetic  field  over 
geologic  time  scales  [Ref.  4] . 
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B.  GEOMAGNETIC  FLUCTUATIONS 

The  geomagnetic  field  changes  with  time.  The  changes  in 
the  field  have  periods  ranging  from  fractions  of  seconds  to 
thousands  of  years.  The  very  slow  variations  (periods  of 
years  or  more)  are  not  of  interest  here  but  are  known  as 
secular  variations  and  are  thought  to  be  geologic  in  origin. 

The  geomagnetic  field  changes  of  interest  here  are  those 
fluctuations  that  occur  over  the  period  of  a  day  or  less. 

It  is  customary  to  view  the  magnetic  field  as  the  stable 
main  field  with  these  comparatively  rapid  fluctuations 
superimposed  on  it.  Figure  2.4  presents  the  spectrum  of 
the  superimposed  fluctuations. 

The  fluctuations  (Figure  2.4)  that  are  superimposed  on 
the  main  field  range  in  intensity  from  .05  to  50  nanoteslas 
and  have  various  causes.  Some  major  contributors  to  these 
fluctuations  are  in  the  ionospheric  currents  overhead  and 
in  magnetic  storms  which  probably  originate  in  solar 
eruptions . 

Finally,  arriving  at  the  point  of  interest  here,  are 
those  disturbances  to  the  earth's  magnetic  field  that 
propagate  through  the  earth's  magnetosphere  as  hydrodynamic 
waves.  These  small  perturbations  cover  a  wide  range  in 
amplitude  and  duration  and  are  known  as  micropulsations. 

Figure  2.5  presents  a  graphical  display  of  these  disturbances. 
Micropulsations  are  thought  to  be  the  result  of  fluctuations 
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100  .01 
Period  (sec) 

Figure  2.5  Field  Strength  of  Micropulsations 
(Campbell,  1966] 

in  the  incident  solar  wind  affecting  wave -particle  interaction 
in  the  magnetosphere  [Ref.  5] .  As  they  propagate  through 
the  conductive  layers  of  the  atmosphere,  they  are  trans¬ 
formed  into  electromagnetic  waves  via  ionospheric  currents, 
which  we  detect  at  the  earth’s  surface  as  magnetic  field 
fluctuations.  There  are  two  major  categories  of  micro¬ 
pulsations:  continuous  and  irregular. 

Smoothly  varying,  periodic  pulsations  are  classified  as 
(Pc)  micropulsations.  Their  distinct  periodicity  is  indica¬ 
tive  of  a  resonance  inside  the  magnetospheric  cavity.  One 
model  proposed  to  explain  this  occurrence  of  smooth,  periodic 
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pulses  is  that  the  resonance  is  a  result  of  an  oscillatory 
hydromagnetic  wave.  This  wave  is  associated  with  groups 
of  ions  which  "bounce"  back  and  forth  in  the  magnetosphere 
[Ref.  6]. 

These  continuous  micropulsations  are  further  classified 
according  to  their  periods  as  follows: 

1.  Pci:  (0.15-5  second  period) 

Known  as  "Pearls,"  these  micropulsations 
are  generated  by  the  cyclotron  instability 
of  the  energetic  protons.  They  have  been  posi¬ 
tively  correlated  with  solar  disturbances  and 
occur  during  daylight  hours  in  the  auroral 
zone  and  during  night  and  early  morning  hours  in 
the  midlatitudes.  Their  average  amplitude  is 
approximately  0.05-0.1  nanotesla  [Ref.  7]. 

2.  Pc2:  (5-10  second  period) 

This  is  a  diurnal  phenomenon  that  shows  some 
positive  correlation  with  solar  activity  and  the 
seasons.  They  usually  decrease  in  their  period 
as  magnetic  activity  increases.  Their  average 
amplitude  is  0.1-1  nanotesla. 

3.  Pc3:  (10-45  second  period) 

This  micropulsation  is  similar  to  Pc2  pulsations 
except  for  the  difference  in  period. 

4.  Pc4:  (45-150  second  period) 

Sunspot  activity  appears  to  have  an  effect  on 
Pc4  pulsations.  The  Pc4  pulsations  also  change 
their  period  seasonally  with  an  average  amplitude 
of  5-10  nanotesla. 

5.  Pc5:  (150-600  second  period) 

These  large-scale  pulsations  typically  occur 
during  morning  and  evening  with  amplitudes  of 
10-100  nanotesla.  Their  duration  shows  a  strong 
geomagnetic  latitude  dependence. 

Irregular  (Pi)  micropulsations  appear  to  be  associated 

with  polar  magnetic  disturbances  and  are  observed  during 

night-time  hours.  The  theory  is  that  these  micropulsations 
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arise  from  disturbances  of  the  ionosphere  due  to  particle 
bombardment  from  the  magnetosphere.  Increased  solar  radia¬ 
tion  levels  also  appear  to  have  an  effect  on  their  existence 
Irregular  micropulsations  are  classified  as  follows: 

1.  Pil:  (1-40  second  period) 

These  pulsations  have  a  small  amplitude  and 
usually  occur  at  night  and  early  morning. 

They  vary  in  intensity  from  0.01-0.1  nanotesla. 

They  demonstrate  a  positive  correlation  with 
auroral  disturbances. 

2.  Pi2:  (40-150  second  period) 

Pi2  pulsation  amplitude  range  is  on  the  order 
of  1-5  nanotesla.  They  usually  occur  during 
the  early  evening  hours  but  may  continue  through¬ 
out  the  night.  The  periods  of  these  pulsations 
decreases  with  increasing  magnetic  activity. 

The  Pc5  and  Pc4  micropulsations  last  several  hours. 
However,  the  Pc3-1  micropulsations  and  irregular  pulsations 
have  a  maximum  duration  of  approximately  one  hour  but  may 
cease  after  only  a  few  minutes  [Ref.  8].  Above  3  Hz,  the 
most  pronounced  field  variations  are  the  earth  ionosphere 
cavity  resonances  of  Schumann  resonances  [Ref.  9] . 

C.  OCEAN  WAVE  EFFECTS 

In  addition  to  all  of  the  geomagnetic  fluctuations 
discussed  above,  if  working  in  the  ocean  environment,  one 
must  take  into  account  hydrodynamically  induced  magnetic 
fluctuations.  These  fluctuations  are  internally  produced 
resulting  from  the  movement  of  a  conductor  (seawater)  in 
an  external  magnetic  field.  The  effects  of  waves  have  been 
calculated  [Ref.  10]  and  are  pictured  in  Figure  2.6. 


Induced  Magnetic  Field  per  Meter  Amplitude 
of  the  Surface  Wave 

[Weaver,  1965] 


Weaver’s  equations  have  been  modified  to  incorporate  bottom 
reflection  and  transmission  with  the  result  that  the  effect 
of  bottom  decreased  the  magnitude  of  the  induced  magnetic 
field  by  less  than  3  dB  for  frequencies  greater  than  0.16  Hz 
in  water  depths  less  than  500  meters  [Ref.  11]. 


III.  COLLECTION  SYSTEM 


A.  DATA  COLLECTION  SYSTEM 

The  data  acquisition  system  is  illustrated  in  Figure  3.1. 
Major  components  consist  of: 

(1)  coil  antenna  sensors  (2) 

(2)  preamplifiers  (2) 

(3)  signal  conditioners  (2) 

(4)  pulse  code  modulation  system  (1) 

(5)  radio  transmitter  (1) 

(6)  radio  receiver  (1) 

(7)  tape  recorder  (1) 

For  the  work  reported  in  this  paper,  two  systems  consisting 
of  components  (1)  through  (5)  were  deployed.  They  wer? 
similar  in  most  respects  with  the  major  exception  that  the 
ocean  equipment  was  enclosed  in  watertight  glass  spheres. 

A  thorough  description  of  the  data  acquisition  system  is 
given  in  Reference  12. 

1 .  Coil  Antenna  Sensors 

Each  coil  antenna  sensor  consists  of  two  continu¬ 
ously  wound  coils  which  are  constructed  of  about  5460  turns 
of  18-gauge  copper  magnet  wire.  Glass  spheres  are  used  to 
maintain  the  watertight  integrity  of  the  coils  and  all  other 
immersed  equipment  for  the  ocean  configuration.  See  Figure 
3.2  for  the  ocean  system.  The  coils  are  mounted  orthogonally 
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instrumentation  subsystem 


Figure  3.1  Data  Collection  System 


.COIL  MAGNETOMETERS 


in  a  nonmagnetic  frame  with  each  coil  having  a  resistance 
of  120  ohms  and  an  inductance  of  approximately  9.31  Henries. 


2.  Preamplifiers 

The  preamplifiers  are  model  13-10A  low-noise  ELF 
amplifiers  manufactured  by  Dr.  Alan  Philips.  The  final 
stage  of  the  amplifier  contains  an  active  low-pass  filter 
which  provides  a  sharp  cutoff  for  frequencies  above  20  Hz. 
Frequencies  above  10  Hz  are  not  examined  in  this  paper.  The 
overall  amplifier  gain  for  inputs  of  less  than  2.5  milli¬ 
volts  is  60  dB  in  power.  The  purpose  of  the  filter  is  to 
remove  60  Hz  noise  and  to  prefilter  the  data  before 
digitization. 

3.  Signal  Conditioners 

The  signal  conditioners  receive  analog  signals  from 
the  preamplifiers,  amplify  them  by  about  30  dB  and  limit 
signals  with  peak  amplitudes  of  7.5  volts  (modifiable)  from 
entering  the  pulse  code  modulation  system. 

4 .  Pulse  Code  Modulation  System 

The  pulse  code  modulation  system  was  designed  and 
manufactured  by  Dr.  Robert  Lowe,  Lowecom,  Inc.  The  system 
offers  a  15-channel  analog  input  capability  with  selectable 
sampling  rates  of  2,  4,  8,  16,  32,  64  and  128  samples  per 
second.  The  PCM  system  uses  a  crystal  oscillator  and  inte¬ 
grated  circuitry  to  develop  the  clocking  pulses,  a  16-channel 
analog  multiplexer,  a  16-channel,  12-bit  analog- to-digital 
converter  as  the  basis  for  the  pulse  coding.  The  crystal 


clock  oscillator  loses  1  bit  in  1  million.  The  clock  pulse 
gates  the  analog  multiplexer,  the  analog- to-digital  converter 
and  additional  circuitry  that  form  the  pulse  code  words.  The 
PCM  output  is  organized  in  frames  with  each  frame  headed  by 
a  synchronization  word  which  is  followed  by  the  pulse -coded 
samples  from  PCM  channels  1  through  15. 

5 .  Transmission  and  Recording 

After  encoding,  the  data  is  transmitted  via  VHF 
radio  to  a  Naval  Postgraduate  School  located  receiver  and 
recorded  on  an  instrumentation  tape  recorder.  For  the  data 
of  this  paper,  a  sample  rate  of  32  samples  per  second  was  used 

B.  DATA  ANALYSIS  EQUIPMENT 

1.  Analysis  Hardware 

The  PCM  data  is  fed  into  a  PCM  decoder  which  places 
the  decoded  PCM  onto  a  9 -track  800  BPI  computer  tape.  This 
computer  tape  is  read  into  the  Naval  Postgraduate  School  IBM 
3033  computer  for  manipulation  and  analysis. 

2 .  Data  Analysis  Software 

The  IBM  3033  is  used  to  convert  the  data  to  usable 
form  and  then  to  perform  any  operations  on  the  data  as  desired 
At  present  the  computer  code  is  written  in  Fortran  IV  and 
converts  the  digital  data  to  voltage,  Fast  Fourier  trans¬ 
forms  the  data,  calculates  power  spectral  densities,  deter¬ 
mines  Stokes  parameters  0  through  3  and  provides  phase  and 
coherence  figures.  Appendix  A  contains  test  verification 
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of  the  code.  Appendix  B  contains  the  code  that  prepares 
data  from  one  site  for  storage  in  the  IBM  3033  Mass  Storage 
System  for  later  recall  and  comparison  with  another  site. 
Appendix  C  contains  the  code  that  prepares  and  manipulates 
the  code  from  a  second  site,  retrieves  from  the  Mass  Storage 
System  that  data  previously  prepared  and  manipulates  the 
two  data  sets  simultaneously. 


IV.  COMPUTER  CODE 


A.  INTRODUCTION 

As  mentioned  earlier,  one  of  the  intentions  of  this 
effort  was  to  develop  a  code  that  would  permit  the  comparison 
and/or  simultaneous  geomagnetic  analysis  of  two  different 
sites.  The  present  data  collection  system  provides  for  the 
recording  of  transmitted  data  from  those  sites  onto  different 
tracks  of  the  instrumentation  tape  recorder.  Since  only  one 
track  of  that  analog  tape  can  be  placed  on  one  digital 
computer  tape,  a  means  of  storing  one  channel  had  to  be 
devised.  As  can  be  seen  by  a  quick  review  of  Appendices  B 
and  C,  the  present  technique  that  allows  this  simultaneous 
manipulation  of  data  from  two  sites  is  to  read  the  computer 
tape  from  one  site  and  to  process  that  data  to  the  point 
where  it  is  ready  for  analysis.  These  results  are  then  stored 
in  the  IBM  3033.  Subsequently,  the  same  routine  (with  modifi¬ 
cations  to  fit  the  parameters  of  the  second  site)  are  applied 
to  data  read  from  the  second  site  computer  tape.  Now  both 
sets  of  data  are  in  the  same  form  and  ready  for  manipulation. 

B.  DATA  ANALYSIS  SOFTWARE 

The  two  programs  (one  for  each  of  two  sites)  are 
essentially  the  same  with  the  exception  that  one  program, 
program  B,  is  a  short  version  of  the  other,  program  A.  Thus, 
the  next  paragraphs  will  discuss  in  detail  the  longer  program. 


This  description  fits  equally  well  to  program  B  except  that 
program  B  ends  after  applying  the  system  transfer  functions 
and  placing  these  results  in  the  Mass  Storage  System. 

1 .  The  Main  Program 


The  main  program  is  divided  into  sections  which 
perform  the  following  functions: 

*  Data  input 

*  Fourier  analysis  of  time  series  data 

*  Application  of  system  transfer  function.  Note:  At 
this  point,  program  B  reads  the  data  into  the  Mass 
Storage  System  and  ends.  Program  A,  at  this  point, 
retrieves  that  same  data  from  Mass  Storage. 

*  Data  averaging 

*  Calculation  and  plotting  of: 

(1)  power  spectral  densities  of  the  individual 
site  geomagnetic  orthogonal  components. 

(2)  orthogonal  component  coherence  of  the  geo¬ 
magnetic  field  of  one  site. 

(3)  Stokes  parameters  0  through  3  of  the  individual 
site  orthogonal  components. 

(4)  single  site  right  and  left  circular  polarization 
power  spectral  densities. 

(5)  single  plane  coherence  of  right  and  left 
circular  polarization  between  sites. 

a.  Data  Input 

The  data  is  read  into  the  IBM  3033  via  the  sub¬ 
routine  "read"  which  was  written  by  Tim  Stanton  of  the 
Naval  Postgraduate  School.  The  "read"  subroutine  reads  the 
PCM  data  from  the  computer  tape  and  converts  it  into  an 
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integer  from  0  to  4095.  The  main  program  than  converts 
this  integer  to  voltage.  See  Appendix  B  or  C. 

b.  Fourier  Analysis 

The  next  section  of  the  code  utilizes  a  sub¬ 
routine  developed  at  the  Naval  Postgraduate  School  and 
called  "fourt."  "Fourt"  Fast  Fourier  transforms  the  voltage 
data  from  the  time  domain  to  the  frequency  domain.  See 
Appendix  A  for  validation  of  this  routine. 

c.  Transfer  Functions 

After  the  Fast  Fourier  transform,  the  program 
applies  the  transfer  functions  for  the  system.  Development 
of  the  transfer  functions  is  contained  in  Reference  12. 

The  transfer  functions  are  applied  up  to  10  Hz  as  straight 
line  functions.  The  data  is  now  in  nanoteslas. 

d.  Power  Spectral  Density 

This  quantity  is  computed  by  multiplying  the 
previously  calculated  spectral  density  summation  by  the 
observation  period  and  dividing  by  the  number  of  data 
groups  (typically  20)  that  the  computer  tape  was  segmented 
into,  providing  an  average  power. 

e.  Coherence 

The  next  steps  of  the  program  determined  for  one 
site  the  geomagnetic  orthogonal  components  Stokes  parameters 
0,  1,  2,  and  3,  and  the  phase  of  the  components  and  their 
coherence.  Next,  right  and  left  circular  polarization 
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parameters  are  computed  for  the  individual  sites  followed 
by  calculation  of  the  coherence  of  these  parameters  between 
the  two  sites. 


C .  PARAMETER  THEORY 

1 .  Stokes  Parameters 

A  useful  tool  for  determining  the  state  of  polariza¬ 
tion  of  a  wave  are  the  Stokes  parameters  0,  1,  2,  and  3 
[Ref.  13]. 

Two  linearly  polarized  plane  waves  with  field 
vectors  in  the  directions  of  and  e2  can  be  described  as: 


*  =  ?  E  e  Ci  -  iwt) 

E1  el  E1  e 


E2  *  e2  E2  e 


(i  t-x  -  iuit) 


These  waves  can  be  combined  to  give  a  general  plane  wave 
propagating  in  the  direction  £  =  kn: 

g  (x,t)  *  (ex  Ex  ♦  £2  E2)  e(i  ^’x  ‘  iu,t)  (1) 


If  E1  and  E2  have  the  same  phase,  (1)  represents  a  linearly 
polarized  wave  with  its  polarization  vector  making  an  angle 

0  *  tan'1(E2/E1)  with  E^  and  with  a  magnitude  E  =  >4:*  +  E* 
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If  E^  and  E 2  have  different  phases,  the  wave  (1)  is 
elliptically  polarized. 

In  the  simplest  case,  circular  polarization,  and 
E£  have  the  same  amplitude  but  differ  in  phase  by  90°.  The 
wave  (1)  is  then: 

£  (x,t)  =  E0  Cex  ±  e2)  e  C^*x  '  iujt)  (2) 

Then  the  components  of  (2)  are: 

E  Qc,t)  *  E  cos  (kz  -  ajt) 

A  U 

(3) 

Ey  (x,t)  *  Eq  sin  (kz  -  tot) 

Thus,  at  a  fixed  po^nt  in  space,  the  fields  (3)  are  swept 
around  circularly  at  the  frequency  w.  For  (e^  +  ic2) ,  the 
rotation  is  counterclockwise  if  the  observer  is  facing  the 
oncoming  wave.  This  wave  is  described  as  having  positive 
helicity.  (e^  -  ie2)  ^as  negative  helicity  or  is  rotating 
clockwise  when  looking  into  the  wave. 

The  two  circularly  polarized  waves  (2)  form  an 
equally  acceptable  set  of  basis  fields  for  description  of 
a  general  state  of  polarization  by  introducing  the  complex 
orthogonal  unit  vectors: 


Then  a  general  representation  of  (1)  is: 


I  (I,t)  •  (E+  t  +  E_  e)  e  Ci  k‘x  '  iBt> 


(4) 


where  E+  and  E_  are  complex  amplitudes.  If  E+  and  E  have 
different  amplitudes  but  the  same  phase,  (4)  is  an  elliptically 
polarized  wave  with  principal  axes  of  the  ellipse  in  the  plane 
defined  by  and 

Thus,  the  polarization  content  of  a  plane  wave  is 
known  if  it  can  be  written  in  the  form  (1)  or  (4)  with  known 
coefficients  (E^  E2)  or  (E+,  E_).  Conversely,  how  can  we 
determine  from  observations  on  a  wave  the  state  of  its 
polarization?  The  answer  lies  in  the  use  of  the  four  Stokes 
parameters.  Their  measurement  or  calculation  determines 
the  state  of  polarization  of  a  wave. 

For  a  wave  propagating  in  the  z  direction,  the  scalar 


products : 


are  the  amplitudes  of  radiation  respectively,  with  linear 
polarization  in  the  x  direction,  linear  polarization  in  the 
y  direction,  positive  and  negative  helicity.  The  squares  of 
these  amplitudes  provide  a  measure  of  the  intensity  of  each 
type  of  polarization.  Phase  information  is  obtained  from 
cross  products.  If  we  assume  the  scalar  coefficients  of 
(1)  and  (4)  to  be  a  magnitude  times  a  phase  factor: 


>i5i  e  =  , 


Ej  “  ax  e  ,  E,  ■  a,  e 


2  2 


c  _  _  _iS+  c  _  _  i6- 
E+  =  a_  e  ,  E_  =  a  e 


and  use  the  circular  polarization  basis  (e+,  £  ) ,  we  have: 
s0  ■  Ie+  •  E|2  +  | e _  •  E|2  ■  a2  +  a2 

Sx  =2Re  [  (e+  •  E+)*  (e_  •  E)J  =  2a+  a_  cos(5_  -  <5  +  ) 

S?  -  2Im[  (e  -  E  )*  (e_  •  EJ]  =  2a.  a_  sin(5_  -  5  +  ) 


2.  Coherence 


If  a  plane  wave  of  frequency  v  is  propagating  in  the 
z  direction,  it  has  components: 

Ex(t)  =  -  i2,tvt 

yt)  -  *2(t)ei,,(t)  ■  i2”vt 

that  are  mutually  orthogonal  and  at  right  angles  to  the 
direction  of  propagation  [Ref.  13].  Then  the  complex 
degree  of  coherence  can  be  written  as  : 


xy 


xy_ 


ST 


xx 


yy 


where  J  , . . .  are  the  elements  of  the  coherency  matrix 

XX 


*  * 

EE  E  E 

a.  a, a 

xx  x  y 

1  12 

*  * 
EE  EE 

a  a7ei<f>1'i4>2  a. 

y  x  y  y 

12  2 

The  complex  factor  p  is  the  complex  degree  of  coherence  and 

xy 

is  a  measure  of  the  correlation  between  components  of  the 


field  vector  in  the  x  and  y  directions. 


D.  SOFTWARE  TESTING 


The  basis  set  for  the  code  calculations  is  that  of  a 
circularly  polarized  wave.  Thus  a  circularly  polarized 
wave  was  chosen  as  the  test  input.  The  wave  was  created  by 
generating  a  4  Hz  sine  wave  for  sites  A  and  B  x  coil  input 
and  4  Hz  cosine  wave  for  sites  A  and  B  y  coil  input.  The 
results  of  this  test  are  valid  only  at  the  input  frequency 
as  a  measure  of  the  code  veracity.  One  of  the  characteris¬ 
tics  of  any  FFT  routine  is  that  small  residuals  are  created 
throughout  the  spectrum  of  interest.  These  residuals  are 
primarily  caused  by  the  finite  observation  time.  These 
small -valued  artifacts  will  be  carried  throughout  all  subse¬ 
quent  parameter  calculations  and  will  appear  in  the  test 
plots.  This  inclusion  of  unexpected  values  at  frequencies 
other  than  the  test  input  of  4  Hz  should  be  considered  as 
the  result  of  the  artificiality  of  the  single  value  test 
input.  Input  of  real  world  data  that  covers  a  full  spectrum 
of  values  obscures  the  "FOURT"  routine  propensity  to  create 
these  artifacts. 

1.  Sites  A  and  B  Orthogonal  Component 
Power  Spectral  Density  Test 


input  to  the  x  coil.  It  follows  that  the  average  power 
spectral  density  is: 


PSD 


xa 


Time 

N 


n 

Z 

1 


* 


Sxa  •  Time 


where  time  is  the  observation  time.  Thus,  for  the  test 
inputs  of  sine  4  Hz  into  the  x  coils  at  sites  A  and  B  and 
the  cosine  4  Hz  into  the  y  coils,  we  expect  and  obtain  (see 
Figures  A.l,  A. 2,  A. 11,  A. 12)  a  peak  magnitude  of  .5  at 
4  Hz  at  each  site. 

2 .  Sites  A  and  B  Orthogonal  Component  Coherence  Test 
Coherence  as  defined  above  is : 


xya 


xya 


/S~  /5~~ 

xa  ya 


for  two  orthogonal  components  where  Sx^a  is  defined  by: 


Time  n  * 

Sxya  j  (Xa> 

With  X  and  Y  our  transformed  inputs  and  considering  the  same 
inputs,  we  expect  to  see  a  coherence  of  magnitude  1  at  a 
frequency  of  4  Hz  and  its  harmonics.  Figures  A. 8  and  A. 18 
demonstrate  the  validity  of  the  coherence  calculations  of 


the  code.  Also  the  value  of  1  for  the  test  case  coherence 
at  frequencies  other  than  4  Hz  is  a  manifestation  of  the 
residual  effect  discussed  above. 

3 .  Sites  A  and  B  Stokes  Parameter  1  Test 

Using  the  definitions  developed  in  preceding 
paragraphs,  Stokes  1  is: 


Sla  =  2a+a_sin  (5_  -  fi +)  =  2  [Sx  -  Sy]  =  0 


for  S  =  S  as  is  the  case  for  this  test  input  and  as  is 
x  y 

shown  in  Figures  A. 4  and  A. 13. 

4.  Sites  A  and  B  Stokes  Parameter  2  Test 

Again  using  definitions  developed  in  the  preceding 
paragraphs,  Stokes  parameter  2  is: 


S_a  =  2a+a_sin  (<5_  -  5  +  )  =  0  for  (6_  -  S  +  )  =  0 


Figures  A. 5  and  A. 15  validate  the  code  for  this  calculation 

5 .  Sites  A  and  B  Stokes  Parameter  3  Test 
Stokes  3  is: 


S3a  =  a+  -  a*  =  -C.707)2  -(.707)2  =  -1 


Figures  A. 6  and  A. 16  validate  the  code  for  this  calculation 
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Stokes  0  is  : 


S  =  a?  +  a2  =  2  [S  a  =  S  . ]  =  2 

oa  +  -  L  xa  xbJ 

(See  Figures  A. 3  and  A. 13.) 

7 .  Sites  A  and  B  Phase  Difference  of  x  and  y  Test 
The  phase  difference  of  sine  4  Hz  versus  cosine  4  Hz 

is  90°  or  1.57  radians.  (See  Figures  A. 7  and  A. 17.) 

8.  Sites  A  and  B  Two  Orthogonal  Component  Power 
Spectral  Density  with  Right  and  Left  Circular 
Polarization  Basis  Set  Test 

The  development  of  the  spectral  densities  for  this 
case  is  similar  to  the  development  of  spectral  density  for 
a  single  component.  Starting  with  the  cross  spectra  density 
definition  for  right  and  left  circular  polarization 

BPxya  *  (x  +  iY)  (right  polarized) 

BM  =  (X  -  iY)  (left  polarized) 

xy  a 

It  follows  that  the  spectral  density  is: 

SPxya  -  iT  l  fBPxya)(BPxya)*  (risht  Polarized) 

SMw.  "  h  "  ,)*  Cleft  polarized) 


The  power  spectral  densities  will  be: 
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N  i  xya 
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where  the  time  is  the  observation  time.  Examples  of  this 

calculation  with  a  test  input  of  4  Hz  right  circularly 

polarized  wave  are  in  Figures  A. 9,  A. 10,  A. 19  and  A. 20. 

9.  Coherence  of  Sites  A  and  B  Right  and  Left  Circular 
Polarization  Basis  set  Test' 

Right  and  left  circular  polarizations  are: 

SP 

COP  =  - abxr  ■  (right  polarization) 

/3F“  /5MT 
a  b 


COM 


SM 


abxy 


/SM  /SM, 
a  b 


(left  polarization) 


where  SPabXy  and  SMabxy  are  the  cross  spectra  densities  of 
the  two  orthogonal  components  at  the  sites  A  and  B. 


SP 


ahrv 


Time 


n 

2  (BP  )  (BP  ) 

v  YV3 '  V  YVfl J 


(right  polarization) 


Cleft  polarization) 


SM  ,  =  iiSSS.  I  (BM  )  (BM  .  ) 

abxy  N  1  K  xya' 1  xyb-7 

The  test  case  of  sine  4  Hz  as  the  xa  and  xb  inputs  versus 
cosine  4  Hz  as  the  ya  and  yb  inputs  yields  a  coherence  of  0 
for  left  circular  polarization  and  1  for  right  circular 
polarization.  (See  Figures  A. 20  and  A. 21.) 


V.  ANALYSIS 


Preliminary,  "eyeball"  analysis  was  conducted  on  data 
collected  over  a  period  from  1122  local  17  Aug  82  to  0639 
local  18  Aug  82.  The  timeframes  inspected  within  this 
window  were : 


TIME  (local) 

DATE 

1122-1254 

17  Aug 

82 

1301-1408 

17  Aug 

82 

2028-2200 

17  Aug 

82 

0121-0251 

18  Aug 

82 

0259-0431 

18  Aug 

82 

0507-0639 

18  Aug 

82 

Original  intentions  were  to  collect  data  over  a  full  24-hour 
period.  However,  the  batteries  at  the  Monterey  Bay  site  went 
dead  after  20  hours.  Of  the  20  hours  of  collected  data, 
eight  hours  appeared  dominated  by  system  noise  and  were  not 
included  in  the  analysis.  However,  the  remaining  data  of 
approximately  eight  hours  spans  a  24-hour  period  and  provides 
evidence  of  some  obvious  trends  that  lend  themselves  to  the 
cursory  analysis  conducted  here.  Representative  plots  of 
the  calculated  parameters  are  in  Appendix  D.  Technical 
Report  NPS  61-83-005  contains  a  complete  set  of  plots  of  the 
calculations  that  span  the  full  20  hours  of  data  collection. 

This  complete  set  of  plots  was  used  in  the  analysis  that  follows 


A.  POWER  SPECTRAL  DENSITIES 

Both  sites  exhibit  a  marked  decrease  in  power  spectral 
density  with  frequency.  There  is  an  obvious  exception  to 
this  observation  at  4  Hz  at  the  land  site.  The  sharp  peak 
at  4  Hz  is  attributed  to  60  Hz  aliasing  which  is  a  result 
of  the  Fast  Fourier  Transform  (FFT)  subroutine.  Fast 
Fourier  algorithm  construction  is  such  that  data  being 
transformed  may  be  convolved  back  onto  itself  [Ref.  15]. 

This  convolution  is  routinely  avoided  by  using  a  sample 
rate  that  is  twice  the  upper  frequency  of  the  domain  of 
interest.  A  sample  rate  of  32  samples  per  second  was  used 
in  this  experiment  and  should  have  precluded  aliasing  of 
the  frequency  of  60  Hz,  especially  when  combined  with  the 
amplifier  filter  cutoff  of  20  Hz.  However,  the  60  Hz 
component  of  the  magnetic  field  was  measured  and  found  to 
be  4  nanotesla/Hz  which  is  10s  higher  than  magnetic  noise 
data  [Ref.  16].  The  20  Hz  filter  cannot  completely  filter 
out  signals  of  this  magnitude.  Thus  the  4  Hz  signal  is  a 
result  of  aliasing  of  the  strong  ambient  60  Hz  signal. 

B .  STOKES  PARAMETERS 

All  Stokes  parameters  demonstrate  reasonable  agreement. 
That  is,  Stokes  0  is  obviously  a*  +  a*  and  Stokes  3  is 
obviously  a*  -  a *  as  shown  in  Figures  D.3,  D.6,  D.13  and  D.16. 


1.  Ocean  Site 


£ 


« 
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The  ocean  site  Stokes  2  and  3  parameters  exhibit  a 
time  dependence.  That  is,  in  the  early  a.m.  hours  the  large 
magnitude  of  Stokes  2  is  the  noticeable  feature.  This 
magnitude  decreases  at  sunrise  and  through  the  day  until  it 
reasserts  itself  after  sunset.  A  review  of  Stokes  parameters 
2  and  3  over  the  total  observation  time  (see  Technical  Report 
NPS-61-83-005)  indicates  that  this  falloff  is  independent  of 
frequency.  The  increase  in  magnitude  of  Stokes  2  after 
sunset  is  also  independent  of  frequency  with  the  exception 
of  some  very  discrete  spectral  components. 

2 .  Land  Site 

The  land  site  Stokes  2  and  3  parameters  show  a  time 
dependence  with  a  shorter  cycle.  That  is,  the  magnitude  of 
Stokes  2  is  large  near  sunset,  followed  by  its  decrease 
during  the  morning  hours.  The  early  afternoon  hours  reflect 
an  increase  in  Stokes  2  until  after  sunset.  Figures  5.1 
through  5.4  show  the  decreased  magnitude  of  Stokes  parameters 
2  and  3  for  the  morning  hours. 

C.  X  AND  Y  COIL  PHASE  DIFFERENCE 
1.  Ocean  Site 

The  use  of  Stokes  parameters  in  the  calculation  of 
the  phase  difference  dictates  that  changes  in  the  time 
series  averaged  Stokes  parameters  should  be  evident  in  the 
x  and  y  component  phase  difference.  This  is  the  case  as 
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Figure  5.1  Stokes  Parameter  2,  La  Mesa,  1122-1254  Local,  17  Aug.  82 


Figure  5.2  Stokes  Parameter  3,  La  Mesa,  1122-1254  Local,  17  Aug.  82 


Figure  5.3  Stokes  Parameter  2,  La  Mesa,  0507-0639  Local,  18  Aug.  82 


Figure  5.4  Stokes  Parameter  3,  La  Mesa,  0507-0639  Local,  18  Aug.  82 
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shown  in  Figures  D.7  and  D.17.  The  phase  difference  from 
sunset  to  sunrise  shows  a  definite  decrease  with  an  "eye¬ 
ball"  average  of  .2  to  .4  radians.  Daylight  hours  show 
a  difference  ranging  from  .5  to  .3  radians  with  the  greatest 
difference  near  local  noon. 

2 .  Land  Site 

The  land  site  phase  difference  reflects  its  Stokes 
parameters  with  a  difference  of  about  .3  radians  in  the 
higher  frequencies  at  all  hours  and  a  difference  of  .1 
radians  during  the  daylight  hours. 

D.  X  AND  Y  COIL  COHERENCE 

1 .  Ocean  Site 

Time  dependence  is  again  a  strong  factor  with  high 
degrees  of  coherence  in  the  ocean  site  x  and  y  coils  being 
greatest  during  hours  of  darkness  and  diminishing  at  sunrise. 
Figures  5.5  and  5.6  are  examples  of  this  phenomenon.  The 
exception  to  this  observation  is  at  local  noon  when  the 
coherence  of  the  ocean  coils  is  again  very  high  but  diminishes 
before  and  after  the  noon  time  data. 

2 .  Land  Site 

The  land  site  x  and  y  coil  coherence  reflects  the 
previous  parameter  cycle  time  dependence.  There  is  a  high 
degree  of  c r'  erence  in  the  lower  frequencies  at  sunrise 
that  is  followed  by  diminishing  coherence  in  the  higher 
frequencies.  Local  noon  shows  an  improved  coherence  at 
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Figure  5.5  Coherence  of  x  and  y  Coils,  Mtry  Bay,  18  Aug.  82,  0121-0251  Local 


Figure  5.6  Coherence  of  x  and  y  Coils,  Mtry  Bay,  18  Aug.  82,  1301-1408  Local 


these  frequencies  but  is  followed  by  a  degraded  coherence 
after  sunset. 

E.  CIRCULAR  POLARIZATION  POWER  SPECTRAL  DENSITY 
There  is  no  marked  contrast  between  right  and  left 

circular  polarization  power  spectral  densities.  Both 
calculations  provide  similar  results. 

F.  COHERENCE  OF  CIRCULAR  POLARIZATION 

No  preference  for  right  versus  left  circular  polarization 
is  discernable.  The  most  obvious  trend  is  that  of  a  high 
degree  of  coherence  in  the  frequency  range  .05  Hz  to  .4  Hz 
as  shown  in  Figures  5.7  and  5.8.  Generally,  coherence  of 
the  two  sites  decreases  thereafter  with  frequency. 


Figure  5.7  Coherence  of  Right  Circular  Polarization  Monterey  Bay/La  Mesa 
1122-1254  Local,  17  Aug.  82 
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VI .  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

In  view  of  the  limited  data  base  and  the  preliminary 
state  of  the  analysis,  only  tentative  conclusions  are 
possible  at  this  point. 

(1)  In  the  frequency  range  .05  to  .4  Hz,  the  ocean  and 
land  data  exhibit  a  high  degree  of  coherence.  At  higher 
frequencies  the  coherence  rapidly  diminishes . 

(2)  Both  polarized  and  unpolarized  signals  were  observed 
with  a  preponderance  of  high  degrees  of  polarization  occurring 
at  the  lower  frequencies.  The  nature  of  the  polarization 
appears  to  be  a  function  of  both  time  and  frequency. 

B.  RECOMMENDATIONS 

The  first  and  obvious  recommendation  is  to  collect  more 
data.  This  is  already  in  the  planning  stage.  The  second  and 
also  obvious  recommendation  is  to  conduct  a  more  thorough  ana¬ 
lysis  of  the  data  on  hand.  This  also  is  already  being  planned 
More  concrete  recommendations  concern  data  processing 

One  90 -minute  real  time  data  tape  requires  12  computer 
cylinders  of  storage  space.  Thus,  for  the  approximately 
eight  hours  of  data  computed  and  analyzed  here,  72  cylinders 
of  storage  were  used.  This  is  too  much  and  is  mostly  a 
matter  of  the  computer  code  employed.  A  little  work  can 
provide  a  program  that  would  require  less  computer  storage 
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space.  This  will  be  a  real  necessity  as  a  data  base  of 
weeks  or  months  is  compiled  for  analysis.  It  is  recommended 
that  the  present  code  by  "streamlined"  to  use  the  IBM 
3033  as  a  calculation  "scratch  pad."  Then  a  few  lines  of 
"Computer  Job  Control  Language"  can  be  employed  to  store 
the  calculations  on  a  high-density  tape  for  future  use  at 
the  Naval  Postgraduate  School  or  anywhere  in  the  world. 


APPENDIX  A 

TEST  PROGRAM 

//JOEM1MB  JOB  (1592,01 29) # • FISHER  SMC. 1399' ,CLASS=C 
//♦MAIN  ORG=NPGVM1.  159  2P, L INES=  (60) 

//♦FORMAT  PB,DDNAME=PL  OT. S YSVECT  R, DEST= LOCAL 

//  EXEC  FRTXCLGP,PARM.  LKED='LIST, MAP, XREF' , REGION. GO=2700K 

//SYSUT1  DD  UNIT=SYSDA,SPACE=(CYL,  (8,8)  ) 

//SYSOT2  DD  UNIT=SYSDA,SPACE=  (CYL,  (8,8)  ) 

//SYSLIN  DD  SPACE= (6080  ,  (80,30) )  ,DNIT=SYSD A 

//FORT. S YS I N  DD  DSN=MS S. SYS3 . NON IMSL . SOURCE ( FOURT) , DIS?= SHR 

//  DD  ♦ 

C  A  LARGE  NUMBER  OF  ARRAYS  ARE  UTILIZED 

C  WITH  THE  INTENTION  OF  MAINTAINING  EACH 

C  BIT  OF  ORIGINAL  DATA  OR  CALCULATION, 

C  BEADY  FOR  RECALL  AT  ANY  TIME  FOR  FURTHER 

C  COMPUTAITON .  EQUIVALENT  ARRAYS  ARE  EMPLOYED 


C 


EECAUSE  THE  PLOTTING  ROUTINE  *DRAWP'  IS 


NOT  ABLE  TO  HANDLE  COMPLEX  ARRAYS. 

ARRAYS  ' ITB 1 f 'RTB* , • ALA3'  AND  *  TITLE  *  ARE 

USED  TO  GENERATE  THE  YERSATEC  PLOTTER  OUTPUT 

OCEAN  ARRAYS  ARE  INDICATED  VIA  THE  SUFFIX 

•0*,  AND  LAND  ARRAYS  ARE  INDICATED  WITH  *L'. 

COMP  LEX*  8  XXO(3  192)  ,YYO  (3192)  f 
C  CSO  (8192)  ,COO(8  192)  ,  ZXO  (3192)  , 

C  CRO(8192),CTO(8  192|  , CUO  (3  1 9 2)  , 

C  CVO  (8192)  ,CWO(3  192)  ,CPO(3192), 

C  ZYO  (819  2)  , COMOL  (8192)  ,CDL  (8192)  , 

C  CUL  (8192)  ,C7L(8  192)  ,CWL(8192)  , 

C  COPOL(8192)  ,  XXL  (819  2)  ,  YYL  (8192)  ,ZXL(8192)  , 

C  ZYL  (8192)  ,CSL(8  192)  ,BPD(3192), 

C  BMO  (8192)  ,BPL(3  192)  ,BML(0192)  , 

C  SPOL  (8192)  ,SMOL  (819  2)  rSPO  (8192)  , 

C  SPL  (8192)  ,SMO(8  192)  ,CTL  (3192)  , 


C  CPL  (8192)  ,CRL(9  192)  ,SML(3192) 

DIMENSION  TIME(  8192)  ,FRE2  (8  192)  , 

C  WORK  (16 3 84)  ,  FRQ  2  (81  92) 

INT  E  GER*4  IT  B  ( 1  2) /I  2*0/ 

REAL *4  RTB(28)  /28*0.  0/, 

C  CROO  (16  3  84)  ,  CFO  0(16  384)  ,3  000  (16  384)  , 

C  COOO  (16  384)  ,  CTD 0(16384)  ,  ZVOO  (16384)  , 

C  CWOO  (16384)  ,CTLL(16384)  , 

C  CPLL  (16  384)  ,CRLL(  16384)  r  33LL  (16384)  , 

C  CULL  (16384)  ,CVLL(16384)  , 

C  CWLL  (16384)  , COM OLL ( 16 384)  , COPOLL  (16384) 
C  SPOO  (16384)  ,SPLL(16384)  , 

C  SMOO  (16384)  ,SMLL(1S384) 

REAL  ALAB  (22)  /f  PDXO’  ,  '  PDI D*  ,  '  ST  10*  , 

C  *  ST20«  ,  •  ST30* , '  STOO'  ,  '  PHAO  »  , 

C  '  GAM  C’  ,  ' PDXL ' ,  * PDYL',»ST1L* , »ST2L», 
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C  COPO I, CO MOL/278  528* (0.0 , 0 . 0 )  / 

DATA  TIME, FREQ,  FRQ2/24576*  0. 0/ 

IFR  A  ME=8  192 

NR=  3  0 

FNS  =  FLO AT  (NR) 

P=0  .  0 

DO  70  Li  =  1 , NR 

THE  DO  LOOP  ENDING  WITH  STATEMENT  70  ENABLES 
THE  PROGRAM  TO  PROCESS  A  LARGE  AMOUNT  OF  DATA 
EY  REPEATING  THE  PROCESS  IN  BLOCKS. THE  DATA 
POINTS  FROM  EACH  RUN  THROUGH  THE  DO  LOOP  ARE 
ADDED  TOGETHER  AND  EVENTUALLY  AVERAGED  BY  THE 
NUMBER  OF  RUNS  THROUGH  THE  DO  LOOP.  ' NR * 
REPRESENTS  THE  NUMBER  OF  DATA  SEQUENCES  TO 
BE  AVERAGED.  1  SEQUENCE  CURRENTLY  EQUALS 
8192  DATA  POINTS  FOR  EACH  CHANNEL  OR  256 


SECONDS  OP  DATA.  FOR  THIS  TEST  SINE  AND 
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COSINE  WAVES  OF  4 HZ  WIL  BE  USED. 

DO  60  J  J=1 , IFRA  ME 

XXO  (JJ)  =  SIN(P*.  7853981631 

YYO  (  JJ)  =COS(P*.  7853  98  16  3| 

XXL  (JJ)  =  SIN(P*.  7853  9816  3) 

YYL  ( JJ)  =COS  (P*.  78  53  98163) 

P=P+  1 

CONTINOE 

THE  FOLLOWING  SECTION  GENERATES  THE  TIME 

AND  FREQUENCY  ARRAYS  AND  NORMALIZES 

THE  INPUT  PCM  DATA  TO  VOLTAGE  FORM  IN 

PREPARATION  FOR  THE  FAST  FOUPIER  TRANSFORM 

TO  THE  FREQUENCY  DOMAIN. 

N=8  1  92 


FN=  FLOAT  (N) 


DELT  AT=  1./32 


T=FN  *DELTAT 

DELT  AF=1  ./T 
DO  20  J=1  ,N 

TIME  (J)  =DELT  AT*  FLOAT  (J> 

FREQ  (J)  =DELTAF* FLOAT  (J) 

FRQ2  (J)  =  ALOG  1 0  (  FREQ  (J) ) 

20  CONTINOE 

THE  NEXT  FD  OR  STATEMENTS  PERFORM 

AN  FFT  ON  THE  INPOT  TIME  SERIES 

DATA. 


CALL  FOORT(XXO,  N,  1,  -1,0,  WORK) 
CALL  FOURT(YYO,  N,  1f  -1,0, WORK) 
CALL  FOORT(XXL,  N,1, -1,0, WORK) 
CALL  FOORT(YYL,  N,  1, -1,0,  WORK) 

[.* 
r  ■ 
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DO  4  C  K4  =  1 , N 


ZXO  (II)  =  ZXO(II)  +  (XX 0  (II)  *CO  N JG  (XXO  (II)  )  ) 
ZYO(H)  =  ZYO  (II)  ♦  (YYO  (II)  *  CO  NJG  (YYO(II))) 
CSO  (II)  =CSO(II)  ♦  (XXO  (II)  * 00 NJG  (YYO  (II)  )  ) 

ZXL  (II)  =  ZXL  (II)  ♦  (XX L  (II)  *CON JG  (XXL  (II)  )  ) 

ZYL  (II)  3  ZYL  (II)  ♦  (YY  L  (II)  *  CO  NJG  (YYL(II)  )  ) 

CSL  (II)  =CSL(II)  +  (XX  L  (II)  *  CO  NJG  (YYL(II))) 

BPO  (H)  =  XXO  (II)  +((0.#1.)  *  YYO  (II)) 

BMO  (H)  =XXO(II)  -(  (0.  ,1.)  *  YYO  (II)  ) 

BPL  (II)  =XXL(II)  +((0.,1.)  *YYL(II)) 

BML  (II)  =XXL  (II )  -  (  (0.  ,1.)  * Y Y L  (II)  ) 

SPOL  (II)  =SPOL  (I  I)  +  (BPO  (II)  *CONJG  (BPL  (II)  )  ) 
SMOL  (II)  =SMOL(I  I)  +  (  BBO  ( II)  *CONJG(EML  (II)  )  ) 
SPO  (H)  =SPO(II)  +  (BPO  (II)  *CONJG  (BPO  (II)  )  ) 
SPL  (II)  =  SPL  (II)  ♦  (BP  L  (II)  *CONJG  (3PL(II)  )  ) 

smo  (n)  =sao(ii)  +(bmo(ii)  *conjg  (bmo  (in  )) 

SML  (H)  =  SML  (II)  +  (BML  (II)  *CONJG  (BML  (II)  )  ) 
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30  CONTINUE 


70  CONTINUE 

THE  NEXT  LOOP  NORMALIZES  THE  ABOVE 
SPECTRA  AND  CALCULATES  POWER  SPECTRAL 
DENSITIES. 

DO  33  13  =  1, N 
ZXO  (13)  =ZXO  (13)  *T/FNR 
ZYO  (13)  =ZY0(I3»  *T/FNR 
CSO  (13)  =CS0(I3)  *T/FNR 
ZXL  (13)  =ZXL(I3)  *T/FNR 
ZYL  (13)  =ZYL(I3)  *T/FNR 
CSL  (13)  =CSL(I3)  *T/FNR 
SPOL  (13)  =SPOL(I  3)  *T/FNE 
SPO  (13)  =SPO  (13)  *T/F  NR 
SPL  (13)  =SPL(I3)  *T/FNR 


SMO  (13)  =SMO  (13)  *T/FNR 


SML  (  13)  =  SML  ( 13 1  *T/FNR 
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SMOL  (13)  =  SMOL(I  3)  *T/FNR 

CONTINUE 

THE  NEXT  LO  OP  CALCULATES  STOKES 

PARAMETERS  0  THROU5H  3  OF  THE 

INDIVIDUAL  SITE  ORTH030NAL 

COMPONENTS  AND  THE  COHERENCE  OF 

THE  PLANAR  CIRCULAR  POLARIZATION 

PARAMETERS  BETWEEN  SITES. 

DO  44  14=1, N 

CTO  (14)  =  (ZXO  (14)  ♦  ZY  0 (14)  »  *2./T 

CPO  (  14)  =  (ZXO  (14  )  -ZYO  (14)  )  *2.  /  (T*CTO  (14)  ) 

CRO  (14)  =  (4*CS0  (  14)  )  /(CIO  (14)  *T) 

COO  (14)  =CSO  (14)  /  (CS QRT  (ZXO  (14)  )  *CSQRI  (ZYO  (14)  )  ) 
CWO  (14)  =  ATAN2  ( A  IM  AG  (COO  (14)  )  ,  REAL  (COO  (14)  )  ) 


COO  (14)  =CSQRT  (C  00  (14)  *COSJS  (COO  (14)  )  ) 


CTO  (14)  =4.34294  48*CLOG  (CTO  (14)  ) 

COO  (14)  =  4.3429448*CLOG(ZXD  (14)  ) 

CVO  (14)  =  4.3429448*CLOG(ZYQ(I4)  ) 

CTL  (14)  =  (ZXL  (14  )  +  ZYL  (14)  )  *2./T 

CPL  (14)  =  (ZXL  (14)  -ZYL  (14)  )  *2./(T*CTL  (14)  ) 

CBL  (14)  =  (4*CSL  (14)  )  /  (CTL(I4)  *T) 

COL  (14)  =CSL  (14)  /  (CSQRT  (ZXL  (14)  )  *CSQRT  (ZYL  (14)  )  ) 

CTL  (14)  =4.34294  48*CLOG(CTL  (14)  ) 

COL  (14)  =4.34294  48*CLOG  (ZXL  (14)  ) 

CVL  (14)  =4.34  294  48*CL0G  (ZYL  (14)  ) 

CWL  (14)  =ATAN2(A  IMAG  (COL  (14)  )  ,  REAL  (COL  (14)  )  ) 

COL  (14)  =CSQRT  (COL  (14)  *CONJG  (COL  (14)  )  ) 
COPOL(I4)=SPOL(I4)/  (CS3RT  (SPO  (14)  )  *CSQfiT  (SPL(I4)  )  ) 
COPO L (14)  =CSQ3T  (COPOL  (14)  *CONJG  (COPOL  (14)  )  ) 
COMOL(I4)  =SMOL(  14) /(CSQRT  (SMO  (14)  )  *CSQRT  (SHL  (14)  )  ) 
C0M0L(I4)  =CS QR P  (COMOL  (14)  *CONJG  (COMOL  (14)  )  ) 


SPO  (14)  =  4.34  294  48 *CLOG  (SPO  (14)  ) 


SPL  (14)  =4. 34  294  48*C LOG  (SPL  (14)  ) 

SMO  (14)  =  4.3429448*CLOG(SSO  (14)  ) 

SML (14)  =4.34294 48*0 LOG (SSL (14) ) 

44  CONTINUE 

C  VERSATEC  PLOT  OF  CALCULATED  QUANTITIES 

NPT  S  =10  .  /DELTA F  +1  . 

C  *NPTS '  DETERMINES  NUMBER  OF  POINTS 

C  NECESSARY  TO  PLOT  THE  0  TO  10  HERTZ  RANGE. 

ITB  ( 2)  =0 
ITB  (  3)  =2  0 
ITB  (4)  =  10 
ITB  (6)  =1 
ITB  (7)  =1 
ITB  (  12)  =1 
RTB  (  1)  =0 
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RTB  (  3)  =  ALAB  ( 1) 


READ  (5,3000)  TITLE 

CALL  DRAWP(RPrS,FRQ2,CUOO,ITB,RTB) 

ITB  (6)  =1 

RTB  (  3)  s  A  LAB  (2) 

READ  (5,3000)  TITLE 

CALL  DRA  WP  (NPTS  ,FRQ2  ,CVO0 , 1 TB,  RTB) 

ITB  (6)  =  1 

RTB  (3)  =  ALAB  (3) 

READ  (5,3000)  TITLE 

CALL  DRAWP(NPTS , FREQ, CP03, ITB, RTB) 

ITB  (6)  =1 

RTB  (3)  =  A  LAB  (4) 

READ  (5,3000)  TITLE 

CALL  DR  A  WP  (NPTS  ,  FREQ,  CROD  ,  ITB,  RTB) 


ITB  ( 6)  =  1 


RTB  (3)  =  ALA3(5) 

READ  (5,  3000)  TITLE 

CALL  DRAWP(NPTS  ,FREQ,CROO  (2)  ,ITB,RTB> 

ITB  <  6)  =  1 

RTB  (  3)  =  ALAB(6) 

READ  (5 , 3 000)  TITLE 

CALL  DRAWP  (NPTS , FRQ2,CTOO ,  ITB,  RTB) 

ITB  (6)  *1 
RTB  (3)  =A LAB  (7) 

READ  (5,  3000)  TITLE 

CALL  DRAWP (NPTS ,FREQ,CWO0 ,ITB, RTB) 

ITB  ( 6)  *  1 
RTB  (3)  =  ALAB{8) 

READ  (5,3  000)  TITLE 

CALL  DRAWP (NPTS , FREQ, COOD, ITB, RTB) 


5 


ITB  (6)  =1 


i 


j 


t*! 


n 

►s 


HTB  ( 3)  =  ALAB  (9) 

READ  (5 , 3000)  TITLE 

CALL  DRAWP (N  PTS  ,  FRQ2  ,  CULL , ITB, RTB) 

ITB  (6)  =1 

RTB  (  3)  =ALAB  (10) 

READ  (5,  3000)  TITLE 

CALL  DRA  WP  (NPTS  ,FRQ2  ,CVLL  ,  ITB,  RTB) 

ITB  (6)  =  1 

RTB  (3)  =  ALAB  ( 1 1 ) 

READ  (5,3000)  TITLE 

CALL  DRAWP(NPT3 , FREQ, CPLL, ITB, RTB) 

ITB  (6)  =1 

RTB  (3)  =ALAB(12) 

READ  (5,  3000)  TITLE 

CALL  DRA  WP  (NPTS,  FREQ,CRLL,  ITB,  RTB) 
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ITB  ( 6)  =  1 


STB  (3)  =ALAB  (13) 

READ  (5,  3000)  TITLE 

CALL  DRAWP{NPTS  ,FREQ,CRLL  (2)  ,ITB,RTB) 

ITB  (6)  =  1 

RTB  (  3)  =  ALAB  ( 14) 

READ  (5,3000)  TITLE 

CALL  DRAWP(SPIS,FRQ2,CTLL, ITB, RTB) 

ITB  (6)  *1 

RTB  (3)  =ALAB  (15) 

READ  (5,3000)  TITLE 

CALL  DRA  WP  (NPTS  ,  FREQ,  CWLL,  ITB,  RTB) 

ITB  (6)  =  1 

RTB  (3)  =ALAB  (16) 

READ  (5,3000)  TITLE 

CALL  DRA  WP  (NPIS  ,  FREQ, COLL,  ITB,  RTB) 
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ITB  (6)  =1 


RTB  ( 3)  =  ALAB  (17) 

READ  (5,3000)  TITLE 

CALL  DRAWP(NPTS , FRQ2, SPOO,ITB,RTB) 

ITB  (6)  =1 

RTB  ( 3)  =ALAB  (13) 

READ  (5, 3000)  TITLE 

CALL  DRA  WP  (NPTS  , FRQ2  , SPLL  ,  I TB,  RTB) 

ITB  (  6)  =  1 

STB  (3)  =  ALAB  (  19) 

READ  (5,3000)  TITLE 

CALL  DRAWP  (NPTS  , FRQ 2 , SM03 , ITB, RTB) 

ITB  (6)  =1 

RTB  (  3)  =A LAB  (20) 

READ  (5,3000)  TITLE 

CALL  DRAW?  (NPTS  ,FRQ2  ,  SMLL  ,  ITB,  RTB) 
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ITB  (6)  =  1 


RTB  (  3)  =  ALAB  (21) 

READ  <5,3000)  TITLE 

CALL  DRA  WP  (NPTS  ,  FEE Q, COPOLL  , IT  B,  RTB) 

ITB  (6)  =1 

RTB  (  3)  =  ALaB  (22) 

READ  (5,3000)  TITLE 

CALL  DRAWP(NPTS , FREQ, COMDLL, ITB, RTB) 
3000  FORM  AT  (6  A8) 

STOP 

END 

/* 

//GO. SYS IN  DD  * 

TEST  POWER  SPECTRAL  DENSITY  OF  SIN  4HZ 
IN  XXO,30  A  VGS ,  (REF  DB  YS  LOG  FREQ) 

TEST  POWER  SPECTRAL  DENSITY  OF  COS  4HZ 
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IN  YYO,30  A  VGS ,  (REF  DB  VS  LOG  FREQ) 


TEST  STOKES  1/STOKES  3  OF  SIN  4HZ 


IN  XXO  AND  COS  4HZ  IN  YYO,  30  AV3S 


TEST  STOKES  2/STOKES  3  OF  SIN  4HZ 


IN  XXO  AND  COS  4HZ  IN  YYOf30  AVGS 


TEST  STOKES  3/STOKES  3  OF  SIN  4HZ 


IN  XXO  AND  COS  4HZ  IN  YYO, 30  AV3S 


TEST  STOKES  0  OF  SIN  4  HZ  IN  XXO  AND 


COS  4HZ  IN  YYO,30  AV3S  ,  (REF  DB  VS  LOG) 


TEST  PHASE  OF  SIN  4  HZ  IN  XXO 


AND  COS  4HZ  IN  YYO,  33  AVGS 


TEST  COHER  OF  SIN  4 HZ  IN  XXO 


AND  COS  4HZ  IN  YYO, 33  AVGS 


TEST  POWER  SPECTRAL  DENSITY  OF  SIN  4HZ 


IN  XXL, 30  AVGS,  (REF  DB  VS  LOG  FREQ) 


TEST  POWER  SPECTRAL  DENSITY  OF  ZOS  4HZ 
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IN  YYL, 30  AVGS. (REF  DB  VS  LOG  FREQ) 

TEST  STOKES  1/STOKES  3  OF  SIN  4HZ 
IN  XXL  AND  COS  4HZ  IN  YYL, 30  AVGS 
TEST  STOKES  2/STOKES  0  OF  SIN  4HZ 
IN  XXL  AND  COS  4HZ  IN  YYL,  30  AVGS 
TEST  STOKES  3/STOKES  3  OF  SIN  4HZ 
IN  XXL  AND  COS  4HZ  IN  YYL,  30  AVGS 
TEST  STOKES  0  OF  SIN  4  HZ  IN  XXL  AND 
COS  4HZ  IN  YYL,  30  AVGS  ,  (REF  DB  VS  LOG) 

TEST  PHASE  OF  SIN  4 HZ  IN  XXL 

AND  COS  4HZ  IN  YYL, 33  AVGS 
TEST  COHSR  CF  SIN  4 HZ  IN  XXL 

AND  COS  4HZ  IN  YYL, 30  AVGS 

TEST  RT  CIRC  POLARIZATION  PSD  OF  SIN  4 HZ  IN 
XXO  AND  COS  4HZ  IN  YYD  ,  (DB  VS  LOG), 30  AVGS 
TEST  RT  CIRC  POLARIZATION  PSD  OF  SIN  4HZ  IN 


XXL  AND  COS  4HZ  IN  YYL,  (DB  VS  LOG)  ,30  AVGS 

TEST  LEFT  CIRC  POLARIZATION  PSD  OF  SIN  4HZ 

IN  XXO  AND  COS  4HZ  IN  YYO  (DB  VS  LOG) ,30  AVGS 
TEST  LEFT  CIRC  POLARIZATION  PSD  OF  SIN  4HZ 

IN  XXL  AND  COS  4HZ  IN  YYL  (DB  VS  LOG) ,30  AVGS 

TEST  COHER  RT  CIRC  POLARIZATION  OF  SIN  4HZ 

IN  XXO, XXL  AND  COS  4HZ  IN  YYO, YYL,  30  AVGS 

TEST  COHER  LEFT  CIRC  POLARIZATION  OF  SIN  4HZ 
IN  XXO, XXL  AND  COS  4HZ  IN  YYO, YYL,  30  AVGS 


G  FREQUENCY  (HZ) 


87 


88 


Figure  A. 6  Test  Ocean  Stokes  Parameter  Three 


FREQUENCY  (HZ) 


Figure  A. 9  Test  Ocean  Right  Circular  Polarization  Power  Spectral  Density 


Figure  A. 10  Test  Ocean  Left  Circular  Polarization  Power  Spectral  Density 
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LOG  FREQUENCY  (HZ) 


AD-A127  708  COHERENCE  STUDIES  OF  GEOMAGNETIC  FLUCTURTIONS  IN  THE 
FREQUENCY  RANGE  05  TO  10  HZ(U>  NAVAL  POSTGRADUATE 
SCHOOL  MONTEREY  CA  J  T  FISHER  DEC  82 
UNCLASSIFIED  F/G  3/14 


Figure  A. 14  Test  Land  Stokes  Parameter  Two 


est  Land  Stokes  Parameter  Three 
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Figure  A. 17  Test  Land  Coherence  of  x  and  y  Coils 
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Figure  A. 19  Test  Land  Left  Polarization  Power  Spectral  Density 


FREQUENCY  (HZ) 

Figure  A. 20  Test  Ocean  Coherence  of  Right  Circular  Polarization 


Figure  A. 21  Test  Ocean  Coherence  of  Left  Circular  Polarization 


APPENDIX  B 


MASS  STORAGE  PROGRAM 

//JOEB1MB  JCB  (  1592,01  29)  ,*  FISHER  SMC. 1399* ,CLASS=F 
//♦MAIN  0RG=NPGVM1.  159  2P 
//  EXEC  FORTXCLG 

//FORT . S YS I N  DD  DSN=HS S. SYS3. NONIMSL. SOURCE (FOURT) , DISP= SHR 
//  DD  ♦ 

C  THE  ARRAY  ‘IN*  MILL  BE  USED  TO 

C  RECEIVE  THE  DATA  PASSED  FROM 

C  THE  SUBROUTINE  •READ*  AND  THEN 

C  TRANSFERRED  TO  THE  APPROPRIATE 

C  XXX  OR  YYY  ARRAY. 

INTEG ER*2  IN(15) 

COMPLEX*8  XXX  (81  92),  YYY(3192) 

DIMENSION  FREQ  (9  192)  ,TIME(8192)  ,  WORK  (16384) 

DATA  XXX,  YYY/ 15 3  84*(  0.0,0.  0)  / 

DATA  TIME, FREQ/1  6384*0.0/ 


c  the  following  section  reads  the 

c  FIRST  SIX  SECONDS  DF  COMPUTER  TAPE 

C  AND  DISCARDS  THIS  DATA. 

ITL=  1  92 

DO  55  JJ*1,ITL 

CALL  RD(20,IN,20  0#IREC,IRR> 

55  CONTI  NOE 
IF  RAM  E*8  192 

NR  *20 

FNR* FLOAT  (NR) 

DO  70  L1*1,NR 

C  THE  NEXT  LOOP  HEADS  THE  COMPOTER 

C  TAPE  OSINS  THE  PROGRAM  PROVIDED 

C  BT  MR.  IIM  STANTON  OF  NAVAL  POST- 

C  GRADOATE  SCHOOL. 

DO  60  JJ* 1 tIFRAM  E 
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CALL  BD(20#IN,10  00#IREC,IRR) 

XXX(JJ)  =  IN  (2) 

TTY ( J  J)  =IN  (3) 

60  CONTI  NOE 

N=81 9  2 

FN=FLOAT  (N) 

DELTA  T*1 . /32. 

T=FN*  DELT  AT 

DELTA  F=1./T 
DO  2  0  J=  1  ,  N 

TIME  (  J)  =  DELT  AT*  FLO  AT  ( J) 

FHEQ  (  J)  =  DELTAF*F  LOAI  (J) 

THE  NEXT  4  STEPS  CONVERT  THE 

DATA  TO  VOLTAGE  AND  ENSURES  THAT 

NO  ERRONEOUS  DATA  HAS  BEEN  INTRODUCED 


INTO  THE  ARRAYS 


XXX(J)=(XXX(J)-2  045.5)  *10./2045.5 


XXX(J)=REAL  <XXX(  J)  ) 

YYY(J)= (YYY(J)  -  2  045.  5) *10.  /2045. 5 
YYY ( J) =RE  AL  (YYY  ( J)  ) 

20  CONTI  NOE 

CALL  FOURT(XXX,N  ,1,-1, 0,  WORK) 

CALL  FOURT(YYY,N,1, -1,0, WORK) 

DO  40  K4=1,N 
XXX(K4)=XXX(K4|  /  FN 
YY Y  ( K 4)  =YYY<K4)  /FN 
40  CONTINUE 

THE  NEXT  LOOP  APPLIES  THE  SYSTEM 
TRANSFER  FUNCTION  TO  THE  TRANSFORMED 


FREQUENCY  DOMAIN  DATA.  THE  TRANSFER 
FUNCTION  CONVERTS  VOLTS  TO  N&NOTESLAS 


DO  9  L=1,N 


FRQ=  F  FEQ  (L) 

1  IF  (FRQ.LE.15.)  30  TO  2 

XXX (L) =XXX  (L)  /(I  05.5-3. 14* FRQ) 

YYY (L) =YYY  (L)  /(I  8 1 . 3  2-7.  59  8 *FRQ) 

GO  TO  8 

2  IF  (FRQ.LE.  10.)  SO  TO  3 
XXX(L)=XXX(L)  /(5 .958*FRQ-30. 97) 

YYY  (L)=YYY  <I)  / (7 . 1 66 * FR2-3 9 .  99) 

GO  TO  8 

3  IF  (FRQ.LE. 7. 5)  SO  TO  4 

XXX  (L)=XXX  (L)  /(3 .4  92*FR2"6.31) 

YYY  (L)=YYY  (L)  /(4 . 2 52*FRQ- 1  0 . 85) 

GO  TO  8 

4  IF  (FRQ.LE.  5.)  GO  TO  5 

XXX  (L)*XXX  (L)  /(2 ,6311*FRQ*0.  14667) 
YYY  (L)=YYY  (L)  /(3  .  0 1 2* FRQ- 1 . 55) 


GO  TO  8 

5  IF  (FRQ.LE.3.)30  TO  6 

XXX  ( L) =XXX  (L)  /(2  .6311*FRQ*0.  14667) 

7  YYY(L)=YYY  (L)  /(2.702*FR2) 

GO  TO  8 

6  XXX(L)=XXX  (L) /(2 .72*FRQ) 

GO  TO  7 

8  CONTINUE 

9  CONTINUE 

THE  NEXT  WRITE  STATEMENTS  SEND 

THE  CONVERTED  DATA  TO  MSS 

FOR  FUTURE  MANIPULATION  AND  RECALL 

THIS  NEXT  SET  OF  STATEMENTS  ALSO 

PROVIDES  THE  USER  A  DIAGNOSTIC. 

WRITE  (21)  XXX 


WRITE  (6,*)  XXX  (1)  ,XXX(8192* 


WRITE  (21)  YYY 


Pi 


WRITE  <6,*)YYY(1|  , YYY  (81921 
70  CONTINUE 

ENDFILE  21 

STOP 

END 

SUBROUTINE  P  D(I  UN ,10  ,  IRS,  IREC,IRQ) 

C 

C  THIS  PROCEDURE  FURNISHED  BY  MR.  TIM  STANTON, 

C  DEPARTMENT  OF  OCEANOGRAPHY. 

C 

C  READ  DATA  FROM  IUN,  ALLIGN  ,  CHECK  &  RETURN 

C 

C  IUN  =TA PE  NUMBER,  EG  20 

C  10=  INTEGER* 2  ARRAY,  16  LONG, 

C  (VALUES  0-4095  ,  SUBTRACT  2048)  *5/2028.  GIVES  VOLTAGE 
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IRS*  NUMBER  DP  PESINCS  ALLOWED  (ERRORS) 


IREC*  COUNTER  OF  RECORDS  (FRAMES  OF  DATA) 
BLOCK  512  BITS,  32  BITS  =  RECORD 
800  BPI  TAPE  CJNLABLED 
IRQ=  NUMBER  DF  ACTUAL  RESINCS  (ERRORS) 


INTEGER  *  2  10  (  16) ,  IP  ( 1 6) 
DATA  IRR  /0/ 

IF  (TREC.EQ.O)  IS=0 
IER-  0 

FORMAT  (16A2) 

IF  (IS.NE.O)  3D  TO  50 
READ  (IUN,20,END=900)  IP 
IREC=IREC+1 

IS=IS+1 


IF  (IS.  LT.17)  SO  TO  50 


READ  (ION, 20, EH  D=900)  IP 
IS=1 

IREC=  IREC  +  1 

ICH=IMASK  (IP  (IS)  ,3,0) +1 
WRITE  (6,55)  IC  H,  IS  , IUN ,  I REC 
FORMAT  (»  RESYNCING  ICH,  IS  ,ION,IREC  ',418) 

IF  (ICH.NE.1)  SO  TO  40 
DO  100  1=1,16 
10  (I )  =ISHIFT  (IP  (IS)  ,4) 

IC H=  I  MAS K  (IP  (IS)  ,3,0)  +1 
IF  (ICH.EQ.I)  SO  TO  80 
IER=IER+ 1 

WRITE  (6,70)  10  N,  IR  EC  ,I,ICH,IER 

FORMAT  ('UNIT',  13, 1  RECORD’  ,16,  'CHAN  S  DATA  CH  »,2I4 


'ERRORS  *  ,17) 


L* 


CONTINUE 


WRITE  (6,130)  IREC,  IRR 

30  FORMAT  (  *  RESYNC  AT  FRAME *, 16 ,* WITH  TOTAL  ERRORS*,  17) 
IER=  0 

IRQ  =  IRR 

GO  TO  50 

50  CONTINUE 

RETURN 

00  WRITE  (6,910)  I  UN, I  REC 

10  FORMAT  (*1  END  OF  UNIT  *,I3,'  AT  REC  *,I7) 

STOP 

END 

FUNCTION  ISHIFT  (IN,NPLC) 

RETURNS  SHIFTED  VALUE  OF  1*2  WORD  IN 


-VE  LEFT,»VE  RIGHT  SHIFT 


INTEGER  *  2  IN 

IP=IN 

IF  (IP.LT.O)  IP=IP*  65536 
IF  (NPLC.LT.O)  GO  TO  30 
ISHIFT-IP/(2**I  ABS  (  NPLC)  \ 

RETURN 

ISHIFT=IP*(2**I  ABS(NPLC) ) 

IF  (ISHIFT.GT.  6  5535)  ISHIFT  =  MOD  (ISHIFT , 65536) 
RETURN 

END 

FUNCTION  IMASK  (IN,  IBL,  IBR) 

MASK  1*2  WORD  IN  OUTSIDE  BITS  IBL  &  IBR 

INTEGER  *  2  IN,  10 

IO=IN 

IF  (IBR.EQ.O)  3  0  TO  50 


IT=  I  SHIFT  (I  N  ,  I B  R) 


10=1  T 

50  IP=ISHIFT  (ID,IBL-15-IBR) 

10=1  P 

IMAS  K=ISHIFT  (ID  ,  1  5-  IBL) 

RETURN 

END 

/* 

//GO.FT21F001  DD  UNII-3330 V,HSV3P=PUBUAr DISP=  (NEW,CATL3)  , 
//  DSN=MSS.  Si 59 2 . GMD1 1  A f 

//  DCB=(RECFM=VBS,BLSSIZE=4096,LRECL=4092) , 

//  SPACE = (C  TL,  (8 , 8) ) 

//GO.  FT20F0  01  DD  UN  IT=  3400 -4  ,  V0L  =  S  ER=GMDT1A,  DISP=  (OLD, PASS) 
//  LABEL=(1  ,NL,  ,IN)  , 


// 


DCB*(RECFM=FB,LREDL=32, BLKSIZE=5 1 2, DEN =2) 


APPENDIX  C 

'  MAIN  PROGRAM 

//JO EK 1MB  JOB  (1592,01  29), 'FISHER  SMC.  1 399 •  , CLASS=  F 

//♦MAIN  ORG=NPGVM1. 159  2P, LINES* (60) 

//♦FORMAT  P3,DD  NAME  =PL  OT. S YSVECI R, DEST=LOCAL 

//  EXEC  FETXCLGP,PARM.  LKED='LISr, MAP, XREF' , REGION. GO=2700K 

//SYSUT1  DD  UNIT=SYSDA,  SPACE*  (CYL,  (8,8)  ) 

//SYSOT2  DD  UNIT=SYSDA,  SPACE*  (CYL,  (8,8)  ) 

//SYSLIN  DD  SPACE*  (6080,  (80,30)  )  ,ONIT=SYSDA 

//FORT . S YS I N  DD  DSN=MS S. SYS3. NON IMSL . SOURCE (FOURT) , DISP* SHR 

//  DD  ♦ 

C  DETAILED  DESCRIPTIONS  OF  THE 

C  VARIOUS  STEPS  AND  LOOPS  ARE 

C  CONTAINED  IN  APPENDICES  A  AND  B. 

INTE GER*2  IN  (15) 

COMP LEX*8  XXO(3  192)  ,YY0  (3  192)  , 

C  CSO(8192)  ,COO(8192)  ,ZX0(9192)  , 


C  CRO  (8192)  ,CTO(8  192)  , COO  (3192), 

C  CVO  (8192)  ,CWO(8  192)  ,CP0  (3192)  , 

C  ZYO  (8192)  r  CO  MOL  (8192)  ,COL  (8192)  , 

C  COL  (81 92)  ,CVL (8  192)  ,  CWL  (3192)  , 

C  COPOL(8  1 92)  ,  XXX  (819  2)  ,XYI  (8192)  , 

C  XXL  (8192)  ,YYL(8  192)  ,ZXL(3192)  , 

C  ZYL  (8192)  ,CSL(8  192)  ,BP3  (819  2)  , 

C  BMO  (8192)  ,BPL(8  192)  ,BML(9192)  , 

C  SPOL  (8192), S MOL  (8192)  ,SPD  (8192)  , 

C  SPL  (8192)  ,SMO(8  192)  ,CTL(8192)  r 
C  CPL  (8192)  ,CRL(8  192)  ,SML(3192) 
DIMENSION  TIME(8192)  ,  PRE3  (8192)  , 

C  WORK  (16  384)  ,  FR 3  2  (81  92) 

INTE GER*4  ITB(1  2)  /I  2*0/ 

REAL  *4  RPB(2  8)  /28*0.  0/,  CROO  (16384) 
C  CPOO  (16384)  ,COOO(16384)  , 


C  COOO  <16384)  ,CT0 0(16 384)  ,CVOO  (16  384)  , 

C  CWOO  (16384)  ,CTLL(  16384)  , 

C  CPLL  (16334)  ,CRLL(16384)  ,  COLL  (16384)  f 
C  CULL  (16384)  ,CVLL(1S384)  , 

C  CWLL  (16  384)  ,  COM OLL(  16384)  ,  COPOLL  ( 16384)  , 

C  SPOO  (16384)  ,SPLL<15384)  , 

C  SMOO  (16384)  ,SMLL(16384) 

REAL  ALAB  (22)  /*  PDXO*  ,  •  PDYO*  ,  •  ST10*  , 

C  'ST20',  »  ST30',  *  STOO*  ,«PHAO*  , 

C  'GAMC* ,'PDXL',*  PDYL ' ,'SIIL' r 'ST2L'f 
C  '  ST3  L*  ,  '  STOL*  f '  PHAL* , *  3  AML'  , 

C  '  SPLO*  ,  *  SPLL '  f  1  SHOO',»SHLL*,'CPOL','CMOLV 
REAL *8  TITLE  (12) 

EQUIVALENCE  (TITLE  (1)  ,RTB<5)  )  ,  (CR00(1)  , 

C  CEO  (  1)  )  ,  (CPOO  (  1)  ,  CP  0(1)  )  , 

C  (C000(1)  ,COO  (1)  )  ,  (CTOO  (1)  ,CTO  (1)  )  ,  (CU00(1)  ,  COO  (1)  )  , 
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C  (CW00(1)  ,CWO  (1 »  )  ,  (CV00(1>  ,  CVO  ( 1)  )  ,  (CTLL  ( 1 )  ,CTL{1)  ) 

C  (CPL  L(1 )  ,CPL(1>  )  ,  (CRLL  ( 1|  ,CRL(1))  ,  (COLL(1)  ,COL(1)  ) 

C  (CtJLL(l)  ,  CUL  (1))  #  (C  V  LL  ( 1 )  rCVL(1)  )  ,  (CHLL(1)  ,CWL(1)  ) 

C  (CO?CLL(1)  ,  COPD  L  ( 1|  )  ,  (COSOLL(I)  , 

C  COMOL(I)  )  ,  (SPDO  (1)  ,SPO(1)  )  , 

C  (SPL  L ( 1 )  ,  SPL  <1 )  )  ,  (S MOO  ( 1 1  ,SMO  (1)  )  ,  (SMLL(1)  ,SML(1>  ) 
DATA  XXO,YYO,CSO,COO,ZXO,ZYO,CRO, 

C  CTO,  CUO,CVO,CW3  ,XXX,YYY, 

C  XXL, YYL,ZXL,ZYL ,CSL,BPO, BMO ,BPL, 

C  BHL,  SPOL,SMOL,  SPO,SPL,SM3, 

C  SML ,CTL,CRL,COL ,CVL, COPOL , CWL, COPOL, 

C  COHOL/294912*(0.0,0.  0)  / 

DATA  TIMS, FREQ, FRQ2/24575 *0. 0/ 

ITL-192 

DO  55  JJ=  1,  ITL 

CALL  RD(20,IN,200,IREC,I£R) 


55  CONTINUE 


I  FRA  ME®  8  192 

NR=2  0 

FN3  =  FLOAT  (NR) 

THE  FIRST  STATEMENTS  OF  LOOP  70 

RECALLS  DATA  PREVIOUSLY  STORED  IN 

THE  IBM  3033  MASS  STORAGE  SYSTEM. 
DO  70  Li =1, NR 
READ  (21)  XXX 
READ  (21)  YYY 
DO  43  II-1,  IFRA  ME 
XXL  (II)  * XXX  (II) 

YYL  (  H)  *YYY  (II) 

43  CONTINUE 

DO  6  0  JJ=1,IFRAME 

CALL  RD  (20,INr1  000, IRECfIRR) 


XXO  (JJ)  =IN(2) 


YYO  (OJ)  =IN(3) 

60  CONTINUE 
N=81  92 
FN=  FLOAT  (N) 

DELTAT=1./32. 

T=FN  +DELTAT 
DELTAF=1  ,/T 
DO  20  J=1,N 

TIME  (J)  =DELT  AT*  FLOAT  (J) 

FREQ  (J)  =DELTAF*  FLOAT  (J) 

FHQ2  (J)  =  ALOG 10  (  FRE2  ( J) ) 

XXO  ( J)  =  (XXO  ( J)  -  20  45. 5)  *5. /2  045. 5 
XXO  (J)  =REAL  (XXO  (J) ) 

YYO  ( J)  =  (YYO  ( J)  -2045.5)  *5.  /2  045.  5 
YYO  {  J)  *REAL(YYO  (J) ) 
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20  CONTINUE 


CALL  FOURT(XXO,  N,1,  -1,3, WORK) 

CALL  FOURT(YYO,  N,  1,  -1,0,  »ORK) 

DO  40  K4  =  1,  N 

XXO  (K4)  =XXO(K4)  /FN 

YYO  (K4)  =YYO  (K4)  /FN 

40  CONTINUE 
DO  9  L= 1 ,N 
FRQ  =  FREQ  (L) 

1  IF(FRQ.  LE.  15. )  30  TO  2 

XXO  (L)  =XXO  (L) /(-8  1.  36  +17.066*FRQ  -  0 . 5389*FRQ*FR Q) 
YYO(L)  =YYO  (L)  /  (  -63.  435  +1  1. 71  06*FRQ-0. 23279*FRQ*F  RQ) 

GO  TO  8 

2  I F ( F  RQ .  LE.10.)  30  TO  3 

XXO  (L)  =XXO  (L) /( -18.  486  ♦  4.9029*FRQ) 

YYO  (L)  =YYO(L) /( -63.  435  +11.7106*FRQ  -0 . 23  279*FRQ*  FRQ) 
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GO  TC  8 

3  IF(FRQ.  LE.7.5)  30  TO  4 

XXO  (L)=XXO  (L)  /  (  3.  94  7*FRQ-  9 . 368) 

YYO  (L)=YYO  (L)  / (  4.  29  5  +  FRQ-  1  1 . 837) 

GO  TO  8 

4  I F (FBQ.  LE.  5. )  30  TO  5 

XXO  (1)  =XXO  (L)  /{ 2.  8105*FPQ-0. 5657) 

YYO  (L)  =YYO  (L)  /  (  2.  80  0»FRQ-.5  0) 

GO  TO  8 

5  XXO  (I)  =XXO  (L)  /(  2.  723*FRQ) 

YYO  (L)  =YYO  (L)  /  (  2.  71  7*FRQ) 

GO  TO  8 

8  CONTINUE 

9  CONTINUE 

DO  3  C  11  =  1,  N 

ZXO  (II)  =ZXO(II)  (XXO  (II)  *  CO  N  JG  (XXO  (II))) 
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ZYO  (H)  =  ZYO(II)  ♦  (YYO  (II)  *CONJG  (YYO  (II)  )  ) 
/CSO(H)  =CSO(II)  ♦  (XX  0  (II)  *CONJG  (YYO  (II)  )) 
ZXL  (II)  =  ZXL(II)  ♦  (XX  L  (II)  *C0N JG  (XXL  (II)  )  ) 
ZYL  (II)  =ZYL  (II)  ♦  (YYL(II)  *CONJG  (YYL(II)  )  ) 
CSL  (H)  =  CSL  (  II)  ♦  (XX  L  (II)  *  CO  NJG  (YYL(II)  )  ) 
BPO  (II)  =  XXO(II)  ♦((0.#1.)*YYO  (II)  ) 

BMO  (H)  =XXO(II)  -  (  (0.  ,  1.)  *YYO  (II)  ) 

BPL  (II)  =XXL  ( II)  ♦(  (0.  ,  1.)  *  YYL  (II)  ) 

3  ML  (H)  =XXL  (II)  -((0.,1.)  *YYL  (II)  ) 

SPOL  (II)  =  SPOL  (I  I)  ♦  (BPO  (II )  *CONJG  (BPL  (II)  )  ) 
SMOL  (II)  =SMOL (I  I)  ♦  (BMO  (II)  *CONJG(BML  (II)  )  ) 
SPO  (H)  =SPO(II)  +  (BPO  (II)  *CONJG  (BPO  (II)  )  ) 
SPL  (  II)  *  SPL  (II )  ♦  (BP  L  (II)  *CONJG  (BPL  (II)  )) 
SMO  (II)  =  SMO  ( II)  ♦  (BMO  (II)  *CONJG  (BMO  (II)  )  ) 
SML  (H)  *  SML  (II)  ♦  (BML(II)  *CONJG  (BHL  (II)  )  ) 

30  CONTINUE 
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70  CONTINUE 


DO  33  13  =  1, N 
ZXO  (13)  =  ZX0(I3)  *T/F NR 
ZYO  (13)  =ZYO(I3)  *T/FNR 
CSO  (13)  =CSO  (13)  *T/FNR 
ZXL  (13)  =ZXL(I3)  *T/F NR 
ZYL  (13)  =ZYL  (13)  *T/FNR 
CSL  (13)  =CSL  (13)  *T/FNR 
SPOL  (13)  =SPOL(r3)  *r/FNR 
SPO  (13)  =SPO(I3)  *T/FNR 
SPL  (13)  =SPL(I3)  *T/FNR 
SMO  (13)  =SH0(I3)  *T/FNR 
SML  (13)  =  SML(I3)  *T/FNR 
SMOL  (13)  =  SM0L (I  3)  *T/FNR 
33  CONTINUE 

DO  44  14  =  1, N 


126 


CTO  (14)  =  (ZXO  (14)  +ZYO  (14)  \  *2./T 


CPO  (14)  »  (ZXO  (14  )  -ZYO  (14)  )  *2.  /(T*CTO  (14)  ) 

C80  (  14)  =  (4*CSD  (  14) )  /  (CIO  (14)  *T) 

COO  (14)  =CSO  (14)  /  (CSQRT  (ZXO  (14)  )  *CSQRT  (ZYO  (14)  )  ) 
CWO  (14)  =ATAN2  (A IMAG  (COO  (14)  )  ,  REAL  (COO  ( 14)  )  ) 

COO  (  14)  =CSQRT  (COO  (14)  *CO» JG  (COO  (14)  )  ) 

CTO  (14)  =4.34294  48*CL0G(CT0  (14)  ) 

COO  (14)  =4.34  294  43*CLOG(ZXO  (14)  ) 

CVO  (14)  =  4.34  294 48*CLOG  (ZYO  (14)  ) 

CTL  (14)  =  (ZXL  (14)  +ZYL(I4>  \  *2./T 

CPL  (14)  =  (ZXL  (14)  -ZYL(I4)  )  *  2  .  /  (T*CTL  (14 )  ) 

CRL  (14)  =  (4*CSL(  14))  /(dL(I4)  *T) 

COL  (14)  =CSL(I4J  /(CSQRT  (ZXL  (14)  )  *CSQRT  (ZYL  (14)  )  ) 
CTL (14) =4.34  294  48*CLOG (CTL (14)  ) 

COL  (  14)  =4.34  294  48*CLOG(ZKL  (14)  ) 

CVL  (14)  =4.34  29  4  48 *C LOG  (ZYL  (14)  ) 


2 


CWL  (14)  =ATAN2(AIMAG  (COL  (14)  )  ,  REAL  (COL  (14)  )  ) 


COL  (14)  =CSQRT  (COL  (14)  *COSJG  (COL  (14)  )  ) 

COPO  L(I4)  =SPOL(I4)/(CSQRr  (SPO(I4)  )  *CSQRT  (SPL(I4)  )  ) 
COPOL(I4)  =CS QR T  (COPOL  (14)  *CONJG  (COPOL  (14)  )  ) 

COMO  1(14)  =SMOL(I4)/(CS3Rr  (SMO  (14)  )  *CSQRT  (SML(I4)  )  ) 
COMOL(I4)  =CSQRT  (COMOL(I4)  *CONJG  (COMOL  (14)  )  ) 

SPO  (14)  =  4.34  294  48*CL0G(SP0  (14)  ) 

SPL  (14)  =4.34294  48*CLOG(SPL  (14)  ) 

SMO  (14)  =  4.34  294 48*CLOG  (SMO  (14)  ) 

SML  (14)  =4.34  294  48*CLOG  (SM  L  (14)  ) 

44  CONTINtJF 


C-'*i 


hi 


f-M ■.  fc  ."I  ■%  777s  ^  ■■.  .1  ,■.  ^  ^  >  s  *.  ^7^  '»C^  .  * 


ITB  (7)  =1 
ITB  (  12)  =  1 
RTB  (  1)  =0 
STB  (2)  =0 
RTB  (  3)  =  ALAB  (1) 

READ  (5,  3000)  TITLE 

CALL  DRAW?  (NPTS  ,  FRQ2  ,COOD  ,ITB,  RTB) 

ITB  (6)=1 

RTB  (  3)  =  ALAB  ( 2) 

READ  (5, 3000)  TITLE 

CALL  DRA  WP (N  PTS  , FRQ2, C700 , ITB , RTB) 

ITB  (6)  =1 

RTB  (  3)  =ALAB  (3) 

READ  (5,3000)  TITLE 

CALL  DR  A  WP  (NPTS  , FREQ, CPOD  ,  ITB,  RT3) 
ITB  (6)  =  1 


— v  mi* ^  --  ---  -*  -  -  -  * 


129 


RTB  (3)  =ALAB  (4) 

READ  (5,3000)  TITLE 

CALL  DRAWP (NPTS ,FREQ,CROO,ITB, RTB) 

ITB  (6)  =1 

RTB  (  3)  =  ALAB  ( 5) 

READ  (5,3000)  TITLE 

CALL  DRAWP(NPTS  ,FREQ,CROD  (2)  , ITB, RTB) 

ITB  (6)  =1 

RTB  (  3)  =A  LAB  (6) 

READ  (5,  3000)  TITLE 

CALL  DRA  WP  (NPTS  ,  FRQ2  ,CT03  ,  IT B,  RTB) 

ITB  (6)  =  1 
RTB  (3)  =  ALAB  (7) 

READ  (5,3000)  TITLE 

CALL  DRAWP (NPTS , FREQ, CHOD , ITB, RTB) 

ITB  (6)  =1 
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STB  ( 3)  *  ALAB  (8) 


READ  (5 , 3000)  TITLE 

CALL  DRA  WP  (NPTS  ,FREQ  ,COOD  ,  ITB,  RTB) 

ITB  (6)  =1 

RTB  (3)  =  A  LAB  (9) 

READ  (5,  3000)  TITLE 

CALL  DR  A  WP  (N  PIS  ,  FRQ2  ,  CULL  ,  I TB,  RTB) 

ITB  (6)  =1 

RTB  (3)=ALAB(10) 

READ  (5, 3000)  TITLE 

CALL  DRAWP(NPIS, FRQ2,CVIL, ITB, RTB) 

ITB  (6)  si 

RTB  (  3)  =ALAB  (11) 

READ  (5,  3000)  TITLE 

CALL  DRA  WP  (NPTS  ,  FREQ,  CPLL,  ITB,  RTB) 
ITB  (6)  =  1 


3 


RTB  (3)  =  A  LAB  (12) 


READ  (5,  3000)  TITLE 

-CALL  DRAWP(NPTS,FREQ,CRLL,ITB,RTB) 

ITB  (6)  =1 

RTB  (  3)  =  ALAB  (13) 

READ  (5,3000)  TITLE 

CALL  DRAWP(NPTS,  FREQ,CRLL  (2)  , ITB, RTB) 

ITB  (6)  =1 

RTB  (3)  =A LAB  (14) 

READ  (5,  3000)  TITLE 

CALL  DRAW?  (NPIS  ,FRQ2  ,CTLL  ,  ITB,  RTB) 

ITB  (6)  =  1 

RTB  (3)  =ALAB(15) 

READ  (5,3000)  TITLE 

CALL  DRAWP(NPIS ,FREQ,CWLL,ITB,RTB) 

ITB  (6)  =1 
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RTB  ( 3)  =  ALAB  ( 16) 


READ  (5,3000)  TITLE 

CALL  DRA  HP (NPTS ,FRE Q, COLL  , ITB,  RTB) 

ITB  (  6)  =  1 

RTB  (3)  =ALAB  (17) 

READ  (5,  3000)  TITLE 

CALL  DR  A  HP  (NPTS  ,FRQ2  ,SPOD  ,  ITB,  RTB) 

ITB  (6)  =  1 

RTB  (3)=ALAB(18) 

READ  (5,3000)  TITLE 

CALL  DRAHP(NPrS , FRQ2, SPLL , ITB, RTB) 

ITB  (6)  =1 

RTB  (  3)  *ALAB  (19) 

READ  (5,  3000)  TITLE 

CALL  DRA  HP  (NPTS  ,  FRQ2  ,  SMOD  ,  ITB  ,  RTB) 
ITB  (  6)  =  1 
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RTB  (3)  =  ALAB(20) 


3000 


READ  (5,3000)  TITLE 

CALL  DRAWP(NPTS  ,  FRQ2  ,  SMLL  ,  ITB,  RTB) 

ITB  (6)  =1 

RTB  (3)  =  ALAB(21| 

READ  (5,3000)  TITLE 

CALL  DRAWP(NPTS , FREQ, COPOLL, ITB, RTB) 

ITB  (6)  =  1 

RTB  (3)  *ALAB(22) 

READ  (5,3000)  TITLE 

CALL  DR  A  WP  (NPTS  , FREQ  , COMOLL  ,  IT  B,  RTB) 
FORMAT  (6 AS) 

STOP 

END 


SOB  ROUTINE  RD  (I  UN,  10,  IRS,  IREC,IRQ) 


c 


THIS  PROCEDURE  FURNISHED  BY  MR.  TIM  STANTON, 


C  DEPARTMENT  OF  OCEANOGRAPHY. 

C 

C  READ  DATA  FROM  IJN,  ALLIGN  ,  CHECK  S 

C 

C  IUN  =TA  PE  NUMBER,  EG  2D 

C  10=  INTEGER* 2  ARRAY,  16  LONG, 

C  (VALUES  0-4095  ,  SUBTRACT  2048)  *5/2028 .  GIVES 

C  IRS=  NUMBER  OF  RESINCS  ALLOWED  (ERRORS) 

C  IRE C=  COUNTER  OF  RECORDS  (FRAMES  OF  DATA) 

C  BLOCK  512  BITS,  32  BITS  =  RECORD 

C  800  BPI  TAPE  UNLABLED 

C  IRQ=  NUMBER-  OF  ACTUAL  RESINCS  (EPRORS) 

C 

c 

INTEGER  *  2  10  (  16)  ,  IP  (16l 


RETURN 


VOLTAGE 
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DATA  IRR  /O/ 

IF  (IREC.EQ.O)  IS=0 
I  ER=  0 

FORMAT  (16A2) 

IF  (IS.NE.O)  3D  TO  50 
READ  (IUN,20,END=900)  IP 
I  REC=IR  EC+1 

IS=IS+1 

IF  (IS.  IT. 17)  30  TO  50 
READ  (IDN,20,END=900)  IP 
13  =  1 

IREC*IREOl 

ICH=I  MASK  (IP  (15)  ,3,0)  +1 
WlvITE  (6,55)  I3H,IS  ,ION,IREC 
FORMAT  (’  RESYNCING  ICH,  IS  ,  ION  ,  IRE 


IF  (IER.EQ.O)  30  TO  150 
IRR=IRR+1 

IF  (IRR.LT.IRS)  GO  TO  120 
WRITE  (6,110) 

0  FORMAT  ('1  STOPPED  IN  SUB  RD  BECAUSE 
C  OF  IRR.GT.', 16,  *  AT  L 1 1 0 •  ) 

IRQ  =  IRR 

STOP 

0  CONTINUE 

WRITE  (6,130)  IREC,  IRR 

FORMAT  ('RESYNC  AT  FRAME'  ,16 WITH  TOTAL  ERRORS',  17) 

IER=  0 

IRQ-IRR 

GO  TO  50 

CONTINUE 


RETURN 


900  WRITE  (6,910)  IUN,IREC 

910  FORMAT  ( •  1  END  OF  UNIT  ',13,'  AT  REC  •  ,17) 

STOP 

END 

FUNCTION  ISHIFT  (IN  ,  NFLC) 

C  RETURNS  SHIFTED  VALUE  OF  1*2  WORD  IN 

C  -VE  LEFT,  +  VE  RIGHT  SHIFT 

C 

INTEGER  *  2  IN 

IP=  I N 

IF  (IP.LT.O)  IP-IP*  65536 
IF  (NPLC.LT.O)  GO  TO  30 
ISHIFT=IP/(2**I  ABS  (NPLC)  ( 

RETURN 

30  ISHIFT=IP*(2**I  ABS  (NPLC)  ) 

IF  ( ISHIFT. GT.  6  5535)  ISHI FT  =  MOD  (ISHI FT ,  65536) 
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RET  0  BN 


END 

FUNCTION  IMASK  (IN,IBL,IBR) 

MASK  1*2  WORD  IN  OUTSIDE  BITS  IBL  &  I3R 

INTEGER  *  2  IN,  10 

10=1  N 

IF  (IBR.EQ.O)  30  TO  50 
IT=I  SHIFT  (IN  ,IB  R) 

10=  I T 

IP=I  SHIFT  (10  ,IB  L-  15  -IBR) 

10=  IP 

I  MAS  K=ISHIFT  (10 ,  15-1814 


RETURN 


//GO. FT20F0  01  DD  ON II* 3400 -4 , 7DL  =  S ER=GMDT1 B , DISP= (OLD, PASS) 
//  LABEL=(1  ,NL,  ,IN)  , 

//  DCB=(RECFM=FB,LRECL=32,BLKSIZE=512,DEN=2) 

//GO.FT21F001  DD  UN  IT*  3330  V,  MS  V  3  P=PUB4  A,  DISP  =  (OLD  ,  KEEP)  , 

//  DSN=MS5. S1592.GMDT1A, 

//  DCB=(PECFM=VBS,BLKSIZE=4096,LRECL=4092) 

//GO.SYSDUMP  DD  SYSOOI=A 
//GO. SYS IN  DD  * 

X  COIL  PSD,  MIRY  BAY,  1  122- 1254  L DC AL , 1 7AUS8 2 , 

20  AVGS,  3 2  S/SEC, 5VOLT ,  (DB  NT**2  VS  LOG  FREQ) 

Y  COIL  PSD,  MTRY  BA Y,  1  122- 1254  LOCAL, 17AUG82, 

20  AVGS,  3 2 S/S EC, 5VOLT ,  (DB  N T**2  VS  LOG  FREQ) 

STOKES  NO.  ONE,  MTRY  BAY,  1122-1254 
LOCAL,  17AUG82,  20  AVGS,  3  2S/SEC,  5VOLT 
STOKES  NO.  TWO,  MTRY  BAY,  1122-1254 
LOCAL,  17AUG82,  20  AVGS,  3  2S/SEC,  5VOLT 


STOKES  NO.  THREE,  MTRY  BAY,  1122-1254 
LOCAL,  17AOG82,  20  AV3S,  3  2S/SEC,  5 VOLT 
STOKES  NO.  ZERO,  MIRY  BAY,  1122-1254  LOCAL, 
17AUG82,  20  AVGS,  32S/SEC,  5 VOLT ,  (LOG  VS  LOG) 
PHASE  OF  X  AND  Y,  MTRY  BAY  1122-1254 
LOCAL,  17A0G82,  20  AV3S,  3  2S/SE:,  5VOLT 
COHER  OF  X  AND  Y  COILS,  MTRY  BAY,  1122-1254 
LOCAL,  17A0G82,  20  AV3S,  3  2S/SE3,  5VOLT 
X  COIL  PSD,  LA  MESA , 112 2- 12545 LOCAL, 17 ADG8 2, 

20  AVGS,  32  S/SEC, 10VOLT,  (DB  NT**2  VS  LOG  FREQ) 
Y  COIL  PSD,  LA  MESA, 1 1  22-1254  LDCA L, 1 7ADG82  , 

20  AVGS,  32  S/SEC, 10VOLT,  (DB  NT**2  VS  LOG  FREQ) 
STOKES  NO.  ONE,  LA  MESA,  1  122012  54 
LOCAL,  17AUG82,  20  AV3S,  32S/SE3,  10VOLT 
STOKES  NO.  TWO,  LA  MESA,  1  122-1254 
LOCAL,  17A0G82,  20  AV3S,  32S/SEC,  10VOLT 


STOKES  SO.  THREE,  LA  MESA,  1122-  1254 
LOCAL,  17AOG82,  20  AV3S,  32S/SEC,  10VOLT 
STOKES  NO.  ZERO,  LA  MESA,  1122-1254  LOCAL, 
17AUG82,  20  AVGS,  32S/SEC,  10VOLT,  (LOG  VS  LOG) 
PHASE  OP  X  AND  Y,  LA  MESA,  1122-  1254 
LOCAL,  17AUG82,  20  AV3S,  32S/SE0,  10VOLT 
COHER  OF  X  AND  Y  COILS,  LA  MESA,  1  122-1254 
LOCAL,  17A0G82,  20  AV3S,  32S/SE0,  10VOLT 
RT  CIRC  POLARIZATION  P  SD ,MTRY  BAY,  1 1 22-1254 
LOCAL,  17AUG82, 20  AVGS,  32S/ SEC, 5 Y OLT ,  (DB  VS  LOG) 
RT  CIRC  POLARIZATION  PSD,  LA  MESA,  1  122-1254 
LOCAL, 17AOG82, 20  AVGS,  32S/SEC , 1 0 VOLT,  (DB  VS  LOG) 
LEFT  CIRC  POLARIZATION  PSD,  MIRY  B AY,  1  122-1 2 54 
LOCAL,  17A0G  82,  20  AV3S  ,32S/SEC,5VOLT,  (DB  VS  LOG) 
LEFT  CIRC  POLARIZATION  PSD,  LA  MESA,  1122-1254 
LOC AL,17AOG82,  20  AVGS,  32S/SEC,  1 0 VOLT,  (DB  VS  LOG) 


COHER  RT  CIRC  POLARIZATION  MTRY  BAY/LA  MESA, 
1122-1254  LOCAL,  17A0382,  20  AV3S ,  32S/SEC 
COHER  LEFT  CIRC  POLARIZATION  MTS Y  BAY/LA  MESA 
1122-1254  LOCAL,  17AU382,  20  AV3S ,  32S/SEC 


ower  Spectral  Density,  Mtry  Bay, 
82,  1122-1254  Local 


Figure  D.4  Stokes  Parameter  One,  Mtry  Bay 
17  Aug.  82,  1122-1254  Local 


Stokes  Parameter  Two,  Mtry  Bay 
17  Aug.  82,  1122-1254  Local 


Stoker  Parameter  Three,  Mtry  Bay 
17  Aug.  82,  1122-1254  Local 


FREQUENCY  (HZ 


Figure  D.10  Left  Circular  Polarization  Power  Spectral  Density 
Mtry  Bay,  17  Aug.  82,  1122-1254  Local 


Figure  D.12  y  Coil  Power  Spectral  Density,  La  Mesa, 
17  Aug.  82,  1122-1254  Local 


Figure  D.13  Stokes  Parameter  Zero,  La  Mesa, 
17  Aug.  82,  1122-1254  Local 


I 
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Figure  D.15  Stokes  Parameter  Two,  La  Mesa 
17  Aug.  82,  1122-1254  Local 


OS/£S 
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FREQUENCY  (HZ) 

Figure  D.16  Stokes  Parameter  Three,  La  Mesa 
17  Aug.  82,  1122-1254  Local 


Figure  D.17  Phase  of  x  and  y  Coils,  La  Mesa, 
17  Aug.  82,  1122-1254  Local 


17  Aug.  82,  1122-1254  Local 


LOG  FREQUENCY  (HZ) 

Figure  D.19  Right  Circular  Polarization  Power  Spectral 
Density,  La  Mesa,  17  Aug.  82,  1122-1254 


LOG  FREQUENCY  (HZ) 

Figure  D.20  Left  Circular  Polarization  Power  Spectral 
Density,  La  Mesa,  17  Aug.  82,  1122-1254 
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